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EXECUTIVE  SUNNARY 

The  present  development  of  fault  diagnosis  theory  is  in  response  to  the 
proliferation  of  large  circuits  and  systems  whose  complexity  precludes  the 
use  of  manual  methods  of  trouble  shooting.  Although  great  success  has  been 
achieved  in  the  automation  of  fault  analysis  for  digital  systems,  the 
automated  diagnosis  of  analog  circuits  and  systems  remains  a  difficult 
problem. 

The  fault  diagnosis  method  discussed  in  this  volume  of  the  progress 
report  is  one  part  of  a  two  phase  automated  test  philosophy.  The  complete 
diagnosis  process  consists  of  D.C.  testing  (Volume  I  of  the  report)  which 
considers  the  nonlinear  aspects  of  the  system  (e.g.  bias  point  determination 
and  multi  frequency  testing  applied  to  the  system  linearized  about  its 
operating  point.  The  former  is  most  suited  to  isolating  hard  faults 
(catastrophic  changes  in  parameter  values  such  as  shorts  and  opens)  while 
the  latter  is  more  suited  to  the  determination  of  more  subtle  variations  in 
component  behavior.  Although  this  volume  discusses  only  the  multifrequency 
testing,  remember  that  both  testing  methods  must  be  combined  to  make  a 
viable  approach  to  the  analog  fault  analysis  problem. 

The  various  methods  of  fault  analysis  by  "simulation  after  test" 
appearing  in  the  recent  literature  have  essentially  two  distinct 
orientations.  The  philosophy  of  the  first  approach,  is  to  equate 
measurement  data  with  a  nonlinear  function  of  ALL  pertinent  parameter  values 
of  the  system.  This  is  in  contrast  to  the  computationally  simpler  method  of 
limiting  the  number  of  faults  allowed,  so  that  the  resulting  equations  are 
linear.  Our  research  focuses  on  the  former  approach  for  the  following 
reasons.  System  component  models  are  often  generated  by  linearizing  the 
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actual  components  behavior  about  an  operating  point.  This  linearization  of 
an  actual  component  often  means  replacing  it  with  a  network  of  idealized 
components,  each  characterized  by  a  parameter  value  which  depends  on  the 
operating  point.  A  good  example  of  this  is  a  transistor  whose  linear  model 
is  generally  a  dependent  source  interconnected  with  several  resistors  and 
capacitors.  If  such  an  operating  point  is  incorrect  but  is  not  identified 

conclusively  by  the  D.C.  analysis,  it  is  possible  for  all  the  parameters  in 

> 

the  model  which  are  affected  by  that  bias  point  to  be  faulty.  In  this  case 
the  a  priori  assumption  that  the  number  of  faults  is  small  is  not 
appropriate.  Also  a  focus  on  complete  parameter  identification  is  useful  in 
the  role  of  failure  prediction.  Periodic  analysis  which  is  capable  of 
identifying  all  component  values  allows  minor  variations  of  each  component's 
behavior  to  be  tracked  over  time  and  its  failure  predicted. 

With  the  above  focus  in  mind  this  report  introduces  the  CCM  (component 
connection  model)  and  describes  two  approaches  to  multi  frequency  diagnosis: 
the  composite  system  transfer  function  approach  and  the  tableau  approach. 
The  two  approaches  are  evaluated  theoretically  and  numerically  with  the 
tableau  approach  shown  to  be  far  superior.  In  the  context  of  the  tableau 
approach  a  diagnosability  set  is  developed  and  upper  and  lower  bounds  on  the 
number  of  test  frequencies  necessary  to  conduct  a  diagnosis  are  given. 

The  report  also  includes  investigations  of  various  other  aspects  of  the 
tableau  fault  diagnosis  equations.  In  particular,  a  detailed  investigation 
of  the  quadratic  nature  of  the  tableau  equations  is  presented.  By 
exploiting  this  quadratic  structure  superior  numerical  algorithms  for  the 
solution  of  the  diagnosis  equations varise.  Examples  of  the  above  techniques 
given  throughout  the  report,  serve  to  illustrate  the  utility  of  the  tableau 
approach. 
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CHAPTER  1 

I 

THE  COMPONENT  CONNECTION  MODEL 


1.  Introduction 

Underlying  the  fault  diagnosis  investigations  which  follow  in 

subsequent  chapters  is  an  interconnected  system  model  called  the 

ccmponent  connection  model,  abbreviated  CCM  Cl, 2,33.  The  purpose  of 

* 

this  chapter  is  to  introduce  the  notation  and  modelling  philosophies  of 
the  CCM. 

Two  sets  of  equations  characterize  the  CCM:  (i)  component 
characteristic  equations  and  (ii)  interconnection  or  topological 
equations.  Basically  all  system  models  are  equivalent.  The  differences 
lie  in^the  way  the  system  information  is  displayed.  The  philosophy  ** 
behind  the  CCM  is  the  separate  descriptions  of  the  component  dynamics 
and  system  topology.  Such  descriptions  result  in  large  numbers  of 
equations  and  unknowns.  In  matrix  form  the  equations  are  very  sparse 
and  hence  admit  the  use  of  sparse  matrix  storage  and  manipulation 
techniques  C43.  For  large  systems,  storage  requirements  are  often  lower 
than  smaller  non-sparse  composite  system  representations.  In  some 
applications  the  use  of  the  CCM  has  produced  dramatic  improvements  in 
the  use  of  computing  resources.  The  most  notable  example  is  the  MARSYAS 
program  which  used  the  CCM  in  a  Space  Shuttle  engine  model.  The  program 
performed  the  simulation  26  times  faster  than  previous  simulations  used 
by  NASA  C53. 
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2.  CCN  Basics 

Figure  1  shows  a  linear  system  consisting  of  N  components.  The 
components  are  shown  explicitly  as  numbered  blocks.  The  section  labeled 
interconnections  contains  only  the  information  describing  how  the 
components  are  connected  to  each  other  and  to  the  outside  world. 

The  system  input  is  the  vector  u  and  the  system  output  is  the 
vector  y.  Each  component  of  the  system  has  an  input  and  output  vector 
a.  and  b.  respectively  (i=1,2,...,N).  The  component  inputs  and  outputs 
will  often  be  handled  collectively  so  for  convenience  define: 

a  “  col  (a,|  /d^/. > ./9||)  (2.1a) 

b  —  col  (b^  ,b^,  • . .  rbj|)  (2.1b) 

The  vectors  a  and  b  are  called  the  composite  component  input  and 
composite  component  output  vectors.  The  composite  component  output,  b, 
is  related  to  the  composite  component  input,  a,  via  the  component 
characteristics.  The  composite  component  input,  a,  is  related  to  both  b 
and  u  Gthe  system  input)  through  the  interconnection  or  topological 
equations.  Finally,  the  system  output,  y,  depends  on  the  composite 
component  output,  b,  and  the  system  input  u. 

3.  The  Component  Equations 

Since  the  components  are  linear  they  can  be  characterized  by  a 
transfer  function  in  the  frequency  domain.  Let  the  transfer  function 
matrix  for  component  i  be  Z^(s,r^)  where: 

(i)  s  is  the  Laplace  transform  variable 
(ii)  r^  is  a  vector  parameter  which  characterizes  component  i 
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The  object  of  the  fault  diagnosis  procedure  developed  in  this  paper  will 
be  to  identify  (diagnose)  the  value  of  r-  (1*1, 2,. ..,N). 

Let  a^(s)  and  b..(s)  be  the  Laplace  transforms  of  the  component 
input  and  output  vectors.  Then  the  component  characteristic  equation 
for  component  i  is: 

b. (s)  *  Z.(s,r.)a.(s)  (3.1) 

Define: 

r  *  col(r1,r2,...#rN)  (3.2) 

Z(s,r)  *  block-diag|z^(s,r^)/Z2(s/r2),...,Zn(s,rn)J  (3.3) 

The  characteristic  relation  for  all  the  components  can  now  be  written 
compactly  as 

b(s)  *  Z(s,r)a(s)  (3.4) 

where  a(s)  and  b(s)  are  the  Laplace  transforms  of  the  composite 
component  input  and  composite  component  output  vectors.  The  matrix 
Z(s,r)  is  called  the  composite  component  transfer  function  matrix.  It 

is  sparse  due  to  its  block  diagonal  structure. 

Figure  2  illustrates  the  relationship  of  the  composite  component 
transfer  function  to  the  entire  system. 
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4.  The  Connection  Equations 

The  connection  equations  model  the  constraints  which  result  when 
the  components  are  interconnected.  These  constraints  are  linear  and 
algebraic  representing  conservation  laws  such  as  Kirchoff's  current  and 
voltage  laws.  Let  u(s)  be  the  Laplace  transform  of  the  system  input 
vector  and  y(s)  the  Laplace  transform  for  the  system  output  vector. 
From  Figure  2  note  that  the  inputs  to  the  interconnections  are  u(s)  and 
b(s)  and  the  outputs  from  the  interconnections  are  a(s)  and  y(s).  This 
leads  naturally  to  the  following  form  for  the  connection  equations: 

a(s)  =  L^b(s)  +  L^uCs)  (4.1a) 

y(s)  =  L21b(s)  +  L22u(s)  (4. 1b) 

where  the  L. ^  are  matrices  whose  dimensions  conform  to  the  given  vector 
quantities.  The  L^  are  called  the  connection  matrices  and  they  are 
generally  very  sparse.  This  occurs  because  most  circuit  elements  and/or 
subsystem  blocks  have  physical  connections  to  only  a  few  "neighbors". 
Figure  3  shows  the  system  schematic  further  updated  to  portray  the  role 
of  the  connection  matrices. 

A  small  example  described  in  the  next  section  best  summarizes  the 
above  development. 

5.  A  CCM  Example 

Figure  4  is  a  circuit  consisting  of  three  components,  two  resisters 
and  a  capacitor.  The  components  are  shown  separated  from  their 
positions  in  the  circuit  to  emphasize  the  difference  between  the 
connection  equations  and  the  component  equations.  The  circuit  input  is 


j»efj fc<a»<*>s 


u  and  the  output  is  y  *  Cy^y.,]*.  Using  the  usual  resistor,  capacitor 
equations  as  component  models  produces  the  following  composite  component 
transfer  function  matrix  equation: 


b^s)  1  /R^  0 

b2<s)  *  0  1  4C2 

b3(s)  0  0 


a1  Cs) 
a2(s> 
a3(s) 


(5.1) 


Mote  that  this  equation  is  independent  of  the  connections.  It  was 
derived  by  relating  the  component  outputs  and  inputs  through  the 
components  dynamics. 

Next  consider  the  connections.  By  using  KVL  and  KCL  on  the 
circuit,  the  following  equations  result: 


a3<s) 


ol  bi(s>  nr 


a1  (s) 


a2(s)  =  10  -1  b2(s)  +  0  u(s) 

.  ...  LO  1  0  .  ...  0 


lb3<s>l 


(5.2a) 


V  /.n  b«(s) 

yi<s)  p  o  u]  1  pn 

y2(s)  [0  1  Oj  b2<s)  +  |oj  u 

-*  b3(s> 


(5.2b) 


Note  that  the  interconnection  equations  did  not  require  any  knowledge  of 
the  actual  components  dynamics. 

Most  notable  in  the  example  is  the  large  number  of  zero  entries  in 
the  Z(s,r)  and  matrices.  In  fact  the  number  of  non-zero  entries  in 
the  system  matrices  usually  Increases  as  the  order  of  the  number  of 


components,  (ie  linearly)  resulting  In  extremely  sparse  matrices  for 
large  systems. 


6.  Simulation  in  the  Frequency  Domain 

As  will  be  seen  the  objective  of  fault  diagnosis  is  to  Identify  the 
parameters  of  a  system  by  observing  the  systems  response  (outputs)  to  a 
given  stimuli  (system  input).  Essentially  this  process  is  the  inverse 
of  the  system  simulation  problem  wherein  system  parameters  and  inputs 
are  given  with  the  goal  of  numerically  computing  system  outputs.  Since 
the  CCN  leads  to  a  highly  efficient  method  for  system  simulation,  it 
appears  reasonable  to  suppose  that  its  use  in  the  fault  diagnosis 
problem  may  share  some  of  the  same  advantages.  Thus  we  present  a  brief 
description  of  the  frequency  domain  simulation  problem  as  an  aid  to 
understanding  the  fault  diagnosis  discussions  to  follow. 

Often,  frequency  domain  simulations  proceed  by  first  constructing  a 
composite  system  transfer  function  matrix  S(s)  and  then  solving  directly 
for  the  system  output  y(s)  according  to 


y(s)  =  S(s)u(s) 


(6.1) 


It  is  possible  to  compute  the  composite  system  transfer  function  matrix 
explicitly  in  terms  of  the  CCN  matrices  as  follows: 


S(s)  *  L22  ♦  L21(I-Z(s,r)L11)”  Z(s,r)L1 


(6.2) 


where  the  inverse  is  taken  as  an  element  in  the  field  of  rational 
matrices  Cl ,9D.  But  this  is  not  the  most  efficient  way  to  compute  y(s) 
using  the  CCN.  Instead  equations  3.4  and  4.1  are  combined  to  form  a 
sparse  tableau  given  in  equation  6.3 
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1(s,r)  -1*1  [a<s>]  T  0  * 

-I  L|.|  b(s)  "  -L^2U(*> 


y(s)  =  L21b(s>  +  L22u(s)  (6.3b) 

The  solution  of  the  sparse  tableau  equations  of  6.3a  are 
substituted  Into  6.3b  to  obtain  the  system  output.  The  advantage  of 
this  method  for  system  simulation  comes  from  the  fact  that  the  tableau 
can  be  solved  using  sparse  matrix  technique.  The  tableau  matrix 


M(s,r> 


Is  very  sparse  since  all  Its  component  matrices  are  highly  sparse. 
Therefore  only  the  Information  on  the  nonzero  entries  need  be  stored. 
The  solution  proceeds  by  using  an  LU-decomposition  algorithm  which 
factors  M(s,r)  as  L(s,r)U(s,r)  where  L(s,r)  is  lower  triangular  and 
U(s,r)  is  upper  triangular.  In  general  the  matrices  L(s,r)  and  U(s,r) 
are  both  sparse  and  can  be  stored  in  the  same  space  allocated  to  the 
storage  of  M(s,r). 

Solution  of  6.3  at  a  single  frequency  is  most  efficiently  done  by 
solving 


by  forward  substitution  and 


by  back  substitution.  All  operations  are  efficiently  executed  using 


! 


sparse  Matrix  techniques  C4D.  Coeputations  over  a  range  of  frequencies. 


s  *  jw,  otj  £  (ii  £  «£,  can  Also  proceed  efficiently  via  a  continuation 
algorithm  C9D. 

As  an  example  equation  6.4  presents  the  tableau  equations  for  the 
circuit  of  figure  4,  where  and  6j  * 


0 

0 

-1 

0 

0 


0  -1 

0  0 

g3  0 
0  0 
0  1 
-1  0 


0  ff 

-1  0 

0  “I 
-1  0 
0  -1 
1  0 


aj  (s)‘ 

*2  <*> 
8j(s) 

b1(s> 

b2(s) 

bj(s) 


‘  0  ‘ 
0 
0 

-uCs) 

0 

0 


(6.4) 


7.  Summary 

In  the  CCM  context  the  equations  describing  the  components  and  the 
equations  describing  their  interconnections  remain  separate  and 
distinct.  It  is  therefore  possible  to  distinguish  the  effect  of  the 
component  dynamics  from  the  effect  of  their  topology.  This  will  prove 
to  be  an  extremely  desirable  feature  of  the  CCM  when  applied  to  the 
fault  diagrams  problem. 

The  use  of  the  CCM  for  system  simulation  has  resulted  in  fast  and 
numerically  well  conditioned  computations  when  sparse  matrix  techniques 


.  ! 

i 

i 

L  I 
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are  employed  C4,5D.  It  Is  not  unreasonable  to  believe  that  similar 
advantages  are  waiting  to  be  achieved  by  the  use  of  the  CCN  In  the  fault 
diagnosis  problem. 


CHAPTER  2 


MULTI FREQUENCY  FAULT  DIAGNOSIS  (THE  COMPOSITE  TRANSFER  FUNCTION  APPROACH) 

1.  Introduction 

This  chapter  has  two  objectives: 

(i)  Introduce  the  Multif requency  Fault  Diagnosis  problem;  and 
Cii)  Highlight  one  of  the  two  major  methods  for  dealing  with  the 
diagnosis  problem,  namely  the  Composite  Transfer  Function 
Approach. 

The  following  example  illustrates  the  idea  of  mult if requency  fault 
diagnosis  and  motivates  a  formal  definition  of  the  problem: 

Suppose  one  desires  to  determine  the  values  of  R  and  L  in  the 
circuit  in  figure  5,  but  direct  measurements  of  the  components  are 
prohibited. 

Let  the  available  data  be  the  measurement  of  the  voltage  across  the 
series  RL  impedance,  y(s),  and  the  known  input  current,  u(s>.  The 
composite  system  transfer  function  of  the  circuit  is: 

y(s)  *  (R  +  sL)u(s)  Cl .1 ) 

Suppose  <i)  s^  is  a  real  test  frequency;  (ii)  u(s^)  =  1;  and  Ciii)  the 
measured  output  is  y^  =  y(s^);  where  the  superscript  M  means  "measured". 
The  evaluation  of  equation  (1.1)  at  s  =  s^  with  the  measured  output  y^ 
produces 
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y*  «  R  +  s.,l 


(1.2) 


where  R  and  L  are  the  unknowns.  Clearly  the  data  at  Is  Insufficient 
to  determine  the  values  of  R  and  L;  however  performing  a  second 
measurement  at  test  frequency  S£  t  s^  yields  the  matrix  equation 


(1.3) 


Equation  (1.3)  uniquely  determines  R  and  L  provided  s^  and  are 
distinct.  This  simple  example  illustrates  the  use  of  multifrequency 
testing  to  identify  circuit  parameters.  With  this  example  in  mind  the 
general  fault  diagnosis  problem  can  be  formulated  as  follows. 

Given  the  following  information: 

(i)  a  linear  network  appropriately  modeled  containing  N 
components,  each  component  characterized  by  a  parameter,  r., 
1=1,2,. . «,N; 

(ii)  a  set  of  test  inputs;  and 

(iii)  a  set  of  test  outputs 

determine  the  parameters,  r.,  i=1,2,...,N.  As  the  simple  example 

illustrated  measurements  at  a  single  test  frequency  may  be  insufficient 
to  determine  alt  the  parameters.  Thus  it  is  necessary  to  make  a  series 
of  measurements  at  enough  different  test  frequencies  to  provide 
sufficient  data  to  uniquely  specify  all  parameters. 

There  are  two  major  approaches  to  multi  frequency  fault  diagnosis 
for  linear  systems  based  on  the  use  of  the  CCM: 


(i>  The  Composite  Transfer  Function  Approach  developed  by  Sen  and 
Saeks  C63  and 

(ii)  The  Tableau  Approach  developed  by  Decarlo  and  Gordan  C73. 

The  purpose  of  this  report  is  to  continue  the  development  of  the 
multi  frequency  fault  diagnosis  theory  in  the  context  of  the  tableau 
approach.  A  comparison  of  the  two  approaches  will  best  illustrate  the 
motivation  for  this  choice.  For  this  reason  the  next  four  sections 
present  a  summary  of  the  composite  transfer  function  approach  as  well  as 
an  example  of  its  use. 

2.  The  Fault  Diagnosis  Equations  (Composite  Transfer  Function  Approach) 
Suppose  a  linear  system  has  the  transfer  function  matrix 
representation  S(s,r)  where  r  is  the  ci rcuit/system  parameter  vector, 
i.e.  r  =  col(r1,r2,...,rn).  The  diagnosis  of  a  linear  circuit/system 
begins  by  measuring  S(s,r)  at  a  set  of  q  distinct  real  test  frequencies, 
s  =  s.,  i=1,2,...,q.  The  measured  data  has  the  form  C6]: 

col  Js  (s.j,  r ),  S<s2,r ),...,  S(sq,r)J 

The  measurement  in  this  form  is  a  matrix  valued  function  of  the 
parameter  vector,  r.  It  is  much  more  convenient  to  convert  it  to  a 
vector  valued  function  through  the  use  of  the  "vec"  function  C8D  of  a 


matrix  defined  as: 
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vec(A)  = 


(2.1) 


where  (i)  A  Is  an  mxp  matrix 

< i i )  A.  is  the  ith  column  of  A 
(iii)  vec(A)  has  dimensions  mpxl 

The  composite  transfer  function  matrix  SCs^r)  C9D  can  be  expressed 
explicitly  in  terms  of  the  CCM  matrix  as: 

S(s.,r)  =  L22  +  L21(I-Z(s.^r)L11)“1  Z(s.,r)L12  (2.2) 

i  *  1,2,...,q.  Let  M  be  the  vector  of  all  the  measured  data  defined 
precisely  by: 


vecCS(s^,r)3 

vecCS(s2,r)D 

vecCS(sq,r)] 


(2.3) 


The  measurement  M  can  also  be  expressed  in  terms  the  CCM  matrices 
directly  using  equation  2.2  and  the  following  identity  C8,103: 

vecCX  Y  Z3  =  CZt  ft  X3  vecCYD  (2.4) 

where  (1)  X,  Y  and  Z  are  any  appropriately  dimensional  matrices; 
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(ii)  the  superscript  t  denotes  matrix  transpose;  *nd 
(iii)  $  is  the  "Kronecker  product,"  i.e.  given  a  pxk  matrix  A  and 
a  mxq  matrix  B  their  Kronecker  product  is  the  pmxkq  matrix 

nils 


A  W  B  £ 
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(2.5) 


Using  equation  2.2  and  2.4  in  equation  2.3  produces  the  fault  diagnosis 
equations: 


vec  CS(s^,r)D 
vec  CS(s2,r>3 


vec  CL223  ♦  CL*2  ft  L21<I-2(s1,r)L11)”13vee  CZfs^r)) 
vec  CL223  ♦  CC21  3  L21«-Z(s2,r)L11)”13vec  CZ(s2,r>3 


A  f (r) 


vec  CS(sq,r)3 


vec  Cl-22]  +  Ct12 


4i<w<,kr)Lii) 


3vec  CZ(sq,r>3 


(2.6) 

where  M  is  the  measurement  vector  and  f(r)  is  an  explicit  nonlinear 
vector  valued  function  of  r.  The  fault  diagnosis  problem  amounts  to 
solving: 

F(r)  in-  Hr)  =  6  (2.7) 

where  0  denotes  the  "zero”  vector. 


t 
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As  mentioned  earlier  the  use  of  multiple  test  frequencies  is 
necessary  since  for  most  practical  cases  a  measurement  of  S(s,r)  at  one 
test  frequency  will  not  provide  sufficient  information  to  determine  r. 
Notice  also  that  the  composite  transfer  function  approach  requires  the 
computation  of  the  composite  transfer  function  matrix  explicitly  (as  in 
equation  2.2).  For  large  systems  this  is  no  easy  matter.  Finally  note 
that  the  fault  diagnosis  equations  are  nonlinear.  Furthermore,  although 
it  is  not  readily  apparent  from  equation  2.6,  the  polynomial  order  of 
f(r)  increases  in  proportion  to  the  size  of  the  system.  A  small  example 
at  the  end  of  this  chapter  will  illustrate  these  points. 

3.  Measure  of  Solvability 

Two  key  questions  concerning  the  fault  diagnosis  equations  F(r)  =  e 

are: 

i)  Does  a  solution  exist? 

ii)  If  a  solution  exists  to  what  extent  is  it  specified  by  M? 

First  consider  the  question  of  existence.  Recall  that  M  represents 
measurements  made  on  an  actual  system  under  test.  If  the  model  used  to 
derive  the  function  F(r)  is  assumed  to  be  precise  it  should  be  possible 
to  simulate  the  actual  system  mathematically  by  using  y  =  f(r).  Since 
any  behavior  of  the  actual  system  can  be  simulated  by  f(’),  there  must 
be  some  value  r  for  which  f(r  )  -  M  =  F(r)  =  e,  provided  M  has  no 
measurement  error.  The  parameter  value  r*  represents  at  least  one 
possible  solution  to  the  equation  F(r)  =  6.  This  discussion  of  the 
existence  of  a  solution  to  the  fault  diagnosis  equation  points  out 
possible  pitfalls  which  should  be  considered  in  the  actual  computation 
of  a  solution.  In  general  there  will  be  some  error  associated  with  the 
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measurement  M.  Furthermore  the  model  used  to  derive  the  function  f(* ) 
used  for  system  simulation  may  not  be  exact.  It  is  therefore  possible 
for  the  fault  diagnosis  equations  to  be  inconsistent.  For  the 
subsequent  discussion,  however,  the  existence  of  a  solution  to  F(r)  =  8 
is  assumed. 

Now  that  F(r)  s  8  has  a  solution,  to  what  extent  is  that  solution 
specified  by  the  measurement,  M?  Under  ideal  circumstances  r  is 
uniquely  specified  but  because  the  fault  diagnosis  equations  are 
nonlinear  multiple  solution  may  exist.  In  spite  of  this  it  is  still 
possible  for  M  to  completely  specify  the  desired  solution  if  the  region 
of  consideration  is  restricted  so  that  the  desired  root  is  an  isolated 
solution.  The  solution,  r*,  is  an  isolated  solution  if  there  is  a 
neighborhood  of  r*  which  contains  no  other  solution  of  M-FCr)  *  8  C19D. 
By  a  neighborhood  of  r*  we  mean  a  set  of  points  -Cr |  ||r*-r||2  <  e>  for 
some  e  >  0,  where  J  J • |  1 2  denotes  the  usual  Euclidean  norm. 

When  is  the  desired  solution  "isolated"  and  if  it  is  not  isolated 
how  much  information  does  the  measurement,  M,  provide  about  the 
solution?  Saeks  and  Sen  C63  answer  these  questions  by  defining  a 
parameter  called  the  "measure  of  solvability".  The  reasoning  behind  the 
measure  of  solvability  is  based  on  the  inverse  function  theorem  which 
requires  the  following  three  definitions  Cl 93 : 

Definition  2. 

The  mapping  F:D  c  Rn  *  Rm  is  Frechet  differentiable  at  x  c  int(D) 
(interior  of  D)  if  there  is  a  linear  operator  A:  Rn  ♦  Rm  such  that 


(3.1) 


lim  <1  / 1 1  h  1 1_)  |  |F(x+h)  -  F(x)  -  Ah||,  *  0 
h  ♦  0  c 

Note:  The  usual  representation  of  this  derivative  is  the  Jacobian 
matrix.  The  derivative  of  F<*>  will  therefore  be  denoted  by  the  matrix 
JFC). 

Definition 

A  mapping  F:  D  c  Rn  ♦  Rn  is  a  homeomorphism  of  D  onto  F(D)  if  F  is 
one-to-one  on  D  and  F  and  F-^  are  continuous  on  D  and  F CD)  respectively. 
Definition  3^ 

The  mapping  F:  D  C  Rn  -*■  Rn  is  a  local  homeomorphism  at  x  e  int(D) 
if  there  exist  open  neighborhoods  U  and  V  of  x  and  F(x)  respectively, 
such  that  the  restriction  of  F  to  V  is  a  homeomorphism  between  U  and  V. 

Clearly  the  solution,  r*  to  the  equation  F(r>  =  6  will  be  isolated 
if  F(r)  is  a  local  homeomorphism  at  r*.  This  is  determined  using  the 
following  C19D: 

Inverse  Function  Theorem:  Suppose  that  F:  D  c  Rn  ♦  Rn  has  a  Frechet 
derivative,  J p < . > ,  in  a  neighborhood  of  x°  e  int(D)  which  is  continuous 
at  xQ  and  that  Jp(xQ)  is  non-singular.  Then  Fisa  local  homeomorphism 
at  xQ. 

With  the  inverse  function  theorem  in  mind  define  the  measure  of 
solvability  C6D  of  the  fault  diagnosis  equation  F(r)  =  e  by 

«(r  )  =  N  -  rank  CJc(r  >D  (3.2) 

O  F  O 

where  <i)  N  is  the  number  of  unknown  parameters, 

(ii)  F:  0  c  RN  -  RN, 

(iii)  Jp(*)  is  the  NxN  Jacobian  Matrix  of  F,  and 
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<iv)  r  is  a  solution  point, 
o 

In  other  words  «(r  )  is  the  dimension  of  the  nullspace  of  J_(r  ).  The 
o  F  o 

Inverse  Function  Theorem  guarantees  that  when  s(r  )  =  0  (J_(r  )  is  non- 

O  r  o 

singular)  then  r  is  an  isolated  solution.  If  «(r  )  #  0  then  the  fault 
o  o 

diagnosis  equations  must  be  augmented  with  6(rQ)  additional  independent 
constraints  to  guarantee  that  rQ  is  an  isolated  solution  of  the 
augmented  equations. 

Given  the  form  of  the  fault  diagnosis  equations  in  equation  2.6  the 
Jacobian  J^(r)  is  the  following  C6D: 


8  <L21U-Z<»1,r>L11>'1»C(dvee  Zd^rH/dr 
«CI+L11<I-Z(*2,r)L1ir1Z<*2,r)3L12)t  8  a^CWt^DL^j'SjCtdvtc  Zd2,r»/dr 


<(CI*L11a-Z(sq,r)L11)"1ZCsq,r))L12)t  8  (L21«-Z(Sq,r)L11)*1)>C(d»vc  Z(*q,r)J/dr 


(3.3) 

It  would  appear  from  its  definition  that  the  measure  of  solvability 
is  a  local  property,  that  is  it  varies  as  the  parameter  values  of  a 
given  system  vary.  Saeks  and  Sen  C63  show  however  that  when  Z(s,r)  is  a 
rational  function  of  r,  6<r)  is  "almost  constant"  and  can  be 
characterized  by  its  generic  value,  6.  "Almost  constant"  means  that 
6(r)  »  6  (a  constant)  everywhere  in  Rn  except  possibly  for  some  values 
of  r  lying  in  an  algebraic  variety  C6D.  An  algebraic  variety  in  Rn  is 
the  intersection  of  the  zero  sets  of  a  finite  set  of  non  constant 
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polynomials  n  variables  C93.  For  example  in  R^  the  solution  of  y-x^«0 

is  an  algebraic  variety.  Imagine  a  function  f(x,y)  that  is  constant 

every  where  except  along  the  locus  of  points  described  by  y-x^=0.  Such 

a  locus  occupies  "zero  area"  in  the  x-y-plane  and  it  is  therefore 

appropriate  to  consider  f(x,y)  to  be  "almost"  constant.  Consider  Figure 

2 

6  which  displays  the  locus  of  the  curve  y-x  =0  in  the  x-y-plane.  If  the 
intuitive  idea  of  "area"  (width  x  length)  is  used  the  curved  line  of 
figure  6  has  zero  thickness  and  consequently  occupies  zero  area  in  the 
plane. 

4.  Solution  of  the  Fault  Diagnosis  Equations 

Although  many  methods  exist  for  the  solution  of  nonlinear  equations 
the  fact  that  an  explicit  form  of  the  Jacobian  of  F(r)  is  known  suggests 

that  the  Taylor  series  expansion  for  F(r)  might  be  useful  for  finding  an 

N  N 

approximation  to  the  solution  point  of  F(r>  *  e.  Let  F:  OCR  *  R  and 

suppose  r  =  r*  is  a  solution  point  of  F(r)  =  e  and  r  is  some  known 

0 

value  close  to  r*.  Express  F(r*)  as  a  Taylor  series  about  the  point  r  . 

0 

F(r*)  =  F(r  )  +  J_(r  )(r*  -  r  )  +  higher  order  terms  (4.1) 
o  r  o  o 

Since  r*  is  a  solution  point,  F(r*>  =  0.  Ignoring  the  higher  order 
terms  because  of  the  closeness  of  rQ  to  r*  produces  the  following 
approximation: 

8  ■  F(r  )  +  J_(r  )(r*  -  r  )  (4.2) 

O  r  O  0 

This  equation  provides  a  means  to  estimate  the  value  of  r*.  That 
estimate  is  a  natural  choice  for  a  new  expansion  point  which  suggests 
the  following  iteration  procedure: 


Figure  6.  Locus  of  the  curve  y 
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JF(r<k))  (r<k+1>  -  r<k)>  *  -  F(r<k>)  (4.3) 
(k) 

where  r  is  the  "kth"  estimate  determined  in  the  previous  iteration 
and  r*k+^  is  the  current  estimate  to  be  calculated.  This  method  is 
called  the  Newton-Raphson  iteration  scheme  and  it  converges  reasonably 
fast  when  r^J  =  rQ  (the  initial  guess)  is  close  to  r*  [12].  If  rQ  is 
not  sufficiently  close  to  the  solution  point  it  is  possible  for  the 
sequence  of  estimates  to  diverge.  Navid  and  Willson  propose  an 
algorithm  based  on  the  following  theorem  which  alleviates  this 
divergence  problem  Cl  3D  s 


Global  Convergence  Theorem; 

Let  F:  D  C  Rn  *  Rn  be  continuously  differentiable  and  have  a 
nonsingular  derivative  F'(x)  *  Jp(x)  on  the  open  set  D,  and  assume  there 
exists  x*  c  0  for  which  F(x*)  *  9.  Then,  given  any  x°  e  D,  there  exist 
a  sequence  {Xk>  such  that  the  iterates 


xk+1  *  xk  -  X.  J-(xk)"1  F(xk),  k  *  0,1,... 


(4.4) 


converges  to  the  point  x*.  Moreover,  there  exist  m,  depending  on  x°, 
such  that  for  k  >m  we  may  take  x^  =  i;  and  hence  achieve  a  faster  rate 
of  convergence  to  x*. 

The  xk  in  this  theorem  diminishes  the  change  between  estimates  and 
stablizes  the  sequence  of  estimates.  When  the  estimates  are 
sufficiently  close  to  the  solution  then  X^  becomes  unity  and  equation 
4.4  becomes  the  Newton-Raphson  iteration  step. 


5.  Example 


To  illustrate  the  use  of  the  Fault  Diagnosis  equations,  apply  them 
to  the  following  circuit:  The  input  is  u  =  V-n  and  the  output  is 
y  *  Cl,  VQDt.  The  circuit  CCM  equations  are 
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The  transfer  function  is 
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where  r  =  (G^,C2,G2,G^)t.  Notice  that  the  calculation  of  S(s,r)  using 
the  CCM  matrix  required  the  inversion  of  (I-Z(s^r)L^). 

For  each  test  frequency,  s  =  s.,  equating  the  transfer  function 
equation  to  its  measurement  yields  two  independent  equations  in  four 
unknowns.  Hence  two  real  test  frequencies  are  both  necessary  and 
sufficient  to  uniquely  specify  the  four  unknown  parameters.  A 
discussion  of  the  problem  of  test  frequency  selection  will  appear  in  a 
later  chapter  in  the  context  of  the  tableau  approach  to  fault  diagnosis. 

To  solve  for  the  unknown  parameters  choose  two  real  test 
frequencies  s^  =  1  and  S£  =  2.  Suppose  that  the  measured  transfer 
function  values  are: 


(5.5) 


This  produces  the  following  nonlinear  vector  fault  diagnosis  equation: 
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To  start  the  iteration  process  choose  an  initial  guess  for  r. 

Normally  the  nominal  (or  design)  values  for  the  components  will  suffice. 

The  first  choice  for  xn  should  be  unity  unless  it  causes  the  estimates 

to  diverge.  A  straightforward  check  for  this  is  to  test  the  inequality 

I  |M-F(r^n+^)  |  <  |  |M-F(r^n^)  |  ig.  Where  IMIj  indicates  the  usual 

Euclidean  norm.  If  it  is  true  continue  the  iteration  process  but  if 

false  repeat  the  step  with  a  new  xn*  calculating  a  new  r*n+^.  Navid 

and  Willson  Cl 3D  decrease  x  by  the  factor  i  unti l  the  r^n+^ 

n  t 


estimate 
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satisfies  the  inequality  test. 

6.  Disadvantages  of  the  Transfer  Function  Matrix  Approach 

The  example  illustrates  some  shortcomings  of  the  Composite  Transfer 
Function  Approach  which  become  significant  when  it  is  -applied  to  large 
systems.  First  the  method  requires  the  computation  of  the  composite 
transfer  function  matrix.  This  involves  the  inversion  of  (I-Z(s.^r)L^) 
for  i=1,2,...,q.  Next,  if  an  explicit  expression  for  Jf(r)  is  used  in  a 
Newton  Raphson  iteration  scheme  then  (I-Z(s,r)L^)  must  be  inverted 
symbolically.  This  is  extremely  cumbersome  for  even  moderately  sized 
networks.  Finally  since  the  CCM  matrix  (the  L.^  and  Z(s,r>>  appear 
imbedded  within  the  transfer  function  expression  it  is  difficult  to 
analyze  the  effect  that  each  has  on  the  fault  diagnosis  equations.  For 
example  test  point  selection  is  equivalent  to  choosing  L^.  It  is 
desirable  to  choose  to  minimize  6  for  a  given  system  but  this  is  not 
an  easy  matter  in  the  Composite  Transfer  Function  context.  The  next 
chapter  will  introduce  the  Tableau  Approach  to  the  fault  diagnosis 
problem.  Although  the  Tableau  Form  of  the  Fault  Diagnosis  equations  are 
equivalent  to  those  of  the  Composite  Transfer  Approach  they  are 
formulated  in  a  context  which  completely  avoids  the  short  comings  listed 


above. 


34 


CHAPTER  3 

TABLEAU  APPROACH  TO  MULTIFREQUENCY  FAULT  DIAGNOSIS 

1.  Introduction 

The  Tabteau  Approach  to  multifrequency  fault  diagnosis  is 

essentially  the  fault  diagnosis  problem  developed  in  Chapter  2  recast  in 
the  context  of  the  component  connection  model  tableau  introduced  in 

Chapter  1.  The  rationale  for  this  new  approach  is  two  fold.  First,  if 
the  use  of  the  CCM  frequency  domain  tableau  permits  the  development  of 

highly  efficient  algorithms  for  simulating  a  systems  behavior  C9D  it 

should  offer  comparable  advantages  in  the  fault  diagnosis  problem  which 
can  be  thought  of  as  simulation  in  reverse.  In  particular  the  tableau 
approach  to  fault  diagnosis  should  prove  to  be  superior  to  the  composite 
transfer  function  approach  discussed  in  Chapter  2. 

Even  more  important  than  any  computational  advantages  which  the 
tableau  approach  might  offer  is  the  additional  insight  into  the  nature 
of  the  problem  to  be  gained  by  this  different  perspective.  The  transfer 
function  matrix  approach  for  example  has  offered  little  insight  into 
problems  such  as  the  selection  of  test  points  or  choice  of  optimum  test 
frequencies.  Intuitively  these  issues  center  upon  the  individual 
characteristics  of  each  component  in  relation  to  how  the  components  are 
interconnected.  In  the  CCM  Tableau  the  component  characteristic 
equations  and  the  connection  equations  remain  separate.  It  is  therefore 
more  reasonable  to  investigate  the  role  each  plays  in  the  problem  in 
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'  this  context  rather  than  through  a  transfer  function  matrix  where  both 

sets  of  equations  are  intermingled  and  hence  indistinguishable.  The 

i 

purpose  of  this  chapter  is  to  detail  the  development  of  the  tableau 
approach  and  to  highlight  its  usefulness. 

2.  Development  of  the  Fault  Diagnosis  Equations 

The  development  of  the  fault  diagnosis  equations  begin  with  the 
frequency  domain  tableau  which  was  introduced  in  equation  6.3  of  Chapter 
1.  For  reader  convenience  the  tableau  is  repeated  here  in  equation  2.1: 

r(s,r>  “I  “I  ra<s)i  r  0  1 

-I  hlJ  [b(s)J  =  [-L12u<s)J  (2*1a) 

y(s)  =  L21b(s)  +  L22U(S)  (2.1b) 

Recall  that  y(s)  is  the  composite  system  output,  u(s>  is  the  composite 
system  input  and  b(s)  and  a(s)  are  the  composite  component  outputs  and 
inputs  respectively.  Also  recall  that  the  components  are  characterized 
by  the  parameter  vector,  r,  via  the  composite  component  transfer 
function  matrix,  Z(s,r)  and  that  the  l.j  matrices  describe  the  systems 
interconnections. 

Suppose  that  a  test  input,  u(s)  is  applied  to  a  circuit/system 

M 

characterized  by  equation  2.1  and  that  the  output,  y  (s)  is  observed. 
The  superscript  "M"  emphasizes  the  fact  that  this  is  a  test  measurement. 

M 

Apply  the  simulation  procedure  in  reverse  by  putting  y  (s)  into  the 
output  equation  (equation  2.1b)  and  gathering  known  quantities  together 


to  obtain: 
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L21b<s>  =  yM(s)  -  L22u<s>  (2.2) 

The  only  unknown  in  equation  2.2  is  the  composite  component  output 

vector  b(s).  A  closer  look  at  the  properties  of  is  necessary  in 

order  to  determine  the  nature  of  the  solution  to  equation  2.2. 

In  any  practical  situation  the  number  of  system  outputs  is  less 

than  the  number  of  component  outputs.  (Measurements  at  essentially 

every  point  within  a  system  makes  the  fault  diagnosis  problem  trivial 

but  it  is  unfortunately  seldom  practical.)  To  be  consistent  with  this 

practical  consideration,  make  the  following  assumption: 

Assumption  1:  For  the  nyxn^  matrix  l_2<j,  where  ny  =  number  of  outputs  and 

n.  =  dimension  of  b(s),  n.  >  n  . 
b  by 

Furthermore  it  is  reasonable  to  make  a  second  assumption  concerning 

L21: 

Assumption  2:  The  matrix  has  full  row  rank. 

The  reasoning  to  support  assumption  2  is  that  any  dependent  row  of 
corresponds  to  an  output  which  can  always  be  expressed  as  a  linear 
combination  of  the  inputs  and  the  remaining  outputs.  As  a  diagnosis 
measurement  such  an  output  provides  no  additional  information  about  the 
system  under  test  and  therefore  can  be  discarded.  Discarding  all  such 
outputs  means  eliminating  all  dependent  rows  of  L^.  The  resulting  l_2.j 
matrix  must  have  independent  rows. 

To  solve  equation  2.2  note  that  assumption  2  insures  that  has  a 
right  inverse  C1A3  and  therefore  a  solution  must  exist.  Unfortunately 
assumption  1  implies  that  L21  must  have  a  non  trivial  null  space  and 
therefore  there  are  infinitely  many  solutions  to  equation  2.2.  The 
solution  to  equation  2.2  can  be  characterized  in  the  following  way  [73: 
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b(s)  *  LglJ  £y*(s)  -  L^2u <s>J  ♦  I  a^(s)  vk  (2.3) 

~R 

where:  i)  is  any  right  inverse  of 

ii)  p  is  the  dimension  of  the  null  space  of 

iii)  The  set  tv^v^/.x/V^  are  basis  vectors  which  span  the 
null  space  of  L^/ 

iv)  The  set  are  arbitrary  (and  possibly 

complex)  functions  of  s. 

With  the  aid  of  the  following  definitions: 

V  *  Cv^ | V2 1  -  -  - 1 vp3  (2.4) 

a(s)  -  ColCa^  (s),<*2(s),...,<*p(s)D  (2.5) 

bQ(s)  =  l'5  [yM(s)  -  L22u(s)J  (2.6) 

equation  2.3  becomes: 

b(s)  *  bQ(s)  +  Vo(s)  (2.7) 

The  vector  bQ(s)  can  be  thought  of  as  a  particular  solution  to 

equation  2.2.  The  vector  V<*(s)  on  the  other  hand  represents  the  null 
space  component  of  the  general  solution.  In  other  words  V<*(s)  can  be 

viewed  as  the  ambiguity  in  determining  b(s)  given  the  input-output 

M 

measurement  Cu(s),  y  (s)3.  The  need  for  multiple  test  frequencies  is  to 
resolve  such  ambiguity  which  is  intrinsically  present  in  any  formulation 
of  the  problem.  Recall  the  simple  example  illustrated  by  figure  5  at 
the  beginning  of  Chapter  2,  which  could  not  be  solved  by  using  a  single 
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test  frequency.  The  above  development  characterizes  this  "insolubility" 
as  the  null  space  of  the  matrix. 

Since  more  than  one  test  frequency  will  usually  be  required  for 
fault  diagnosis,  let  q  be  the  number  of  test  frequencies  required. 
Denote  the  test  frequencies:  s  =  s.,  i=1,2,...,q.  The  inverted  output 
equation  (equation  2.7)  becomes  a  family  of  equations: 

b(s.)  =  b0(s.)  +  V«(s.)  i  =  1,2,...,q  (2.8) 

Next  substitute  b(s^)  from  equation  2.8  into  the  connection  equations, 

a(s)  *  L^^b(s)  +  L12^(s),  which  is  part  of  the  tableau  in  equation  2.1a. 
This  provides 

a(s^)  »  £b0*si*  +  Va<s-j)J  +  L^u(s)  (2.9a) 

a(s.)  -  L^bQ(s.)  +  L2^u(s.)  +  L^Va(s.)  i  *  1,2,...,q  (2.9b) 

Since  the  first  two  terms  on  the  right  hand  side  of  equation  2.9b  are 
known,  a(s.)  can  also  be  written  as  the  sum  of  a  particular  solution  and 
an  ambiguity  component: 

a(s.)  *  aQ(s.)  +  L^VoKs.)  i  =  1,2,...,q  (2.10) 

where 

aQ(s^)  =  L^b^s.)  +  L21u(s.)  (2.11) 

Now  substitute  equations  2.8  and  2.10  into  the  remaining  part  of  the 
tableau  equation  (component  characteristic  b(s)=Z(s,r)a(s>>. 
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b0(s.)  ♦  VoKs^  =  Z(s.,r)[a0<s.>  +  L^VmCs^j  i  *  1,2,..,q  (2.12) 


The  unknowns  in  equations  2.12  are  r  and  a(s.j),  1=1,2,. ..,q.  Finally 

rewrite  equation  2.12  in  the  following  product  form  C7D: 


[z(s.,r)  |  -vj 


Rivals.)  +  ao(Si)' 


a(s. ) 


=  bQ(s.)  i  =  1,2,. .,q  (2.13) 


The  family  of  equations  defined  by  equation  2.13  constitute  the  tableau 
form  of  the  fault  diagnosis  equations.  Notice  that  this  set  of 
equations  does  not  require  the  computation  of  the  composite  transfer 
function  matrix.  The  nonlinear  component  of  these  equations  are  the 
terms  Z<s.,r)L11Va(s.j ) .  The  order  of  this  nonlinearity  does  not 
increase  as  the  size  of  the  system  increases.  In  fact  for  systems  where 
the  components  have  the  form  Z.(s,r.)  =  r-Z.(s)  the  equations  are 
quadratic  regardless  of  system  size.  Further  discussion  of  these  points 
follows  in  a  later  chapter. 


3.  Solution  of  the  Fault  Diagnosis  Equations 

The  product  form  of  equation  2.13  lends  itself  to  the  following 
compact  characterization  of  the  fault  diagnosis  equations  C7]: 

=  e  (3.1) 

(1)  Of  *  a(s. ) 


F(x) 


f1<r)g1(o1)-B1 


f  (r)g  (a  >-B 

.  q  q  -q  ql 


where 
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< 1 i )  x  =  col(a1,ag,...,a^,r) 
(iii)  f.(r)  =  CZ(s^,r)|-V3 


n-11Vii+ao<si)‘ 


Civ)  g..^)  = 


a . 

—1 


(v)  8.  =  b  (s.) 
i  o  1 

The  objective  is  to  find  the  value  of  x  for  which  FCx)  =  0.  One 
method  is  to  apply  the  Newton-Raphson  scheme  introduced  in  the  previous 
chapter.  In  the  case  of  the  tableau  form  the  iteration  step  is: 


JF(x<k>)(x<k+1)  -  x<k))  *  -  F(x(k>) 


(3.2) 


CIO 

where  (i)  x  is  the  k-th  estimate  of  the  root 
and  (ii)  Jf  is  the  Jacobian  of  F(.) 

As  a  result  of  the  product  structure  of  equation  3.1  the  Jacobian 
has  a  sparse  elegant  structure  C73: 


JF(x) 


39. 

fl<r>157<-1> 


3g2 


3f  . 


3g_  3f 

f  (r) .■;"(«  )  T"^r*0fl(<O 

q  ao^  — q  3r  q  — q 


(3.3) 


where 


41 


8g. 


f1(r)  837  (“1>  =  z<Vr>  L11  V  “  V 


<3. 4a) 


and 


a 7i(r>  W  = 


(3.4b) 


•"l  !^Vr>  [L11  V  if  *  W] 


i 

Note  that  f.(r)-^— <a.)  Is  independent  of  a..  Also  notice  that  when  the 

I  0CU  “I  —1 

components  have  the  form  Z.(s/r.)  =  r.Z.(s)  (e.g.  RLC  circuit  elements) 
8f .  )  J  )  ) 

then -r~(r)g.  (a.)  forms  a  diagonal  block  of  J_(.)  which  is  independent 

o  i  1  "■!  r 

of  r.  Thus  the  Jacobian  associated  with  the  tableau  fault  diagnosis 
equations  has  an  elegant  simple  sparse  structure  which  compares 
favorably  with  the  "ugly"  structure  of  the  Jacobian  in  the  composite 
approach  (equation  3.3  in  Chapter  II). 


4.  An  Example 

The  example  presented  below  is  the  same  example  which  was  presented 
in  Chapter  2  (figure  7).  Besides  illustrating  the  use  of  the  fault 
diagnosis  equations  in  the  CCM  context  this  will  provide  a  basic 
comparison  of  the  two  approaches.  Such  a  comparison  will  illustrate  the 
relative  ease  with  which  the  fault  diagnosis  equations  are  solved  as 
well  as  the  simplicity  in  computing  the  Jacobian  necessary  for  the 
Newton-Raphson  iteration.  Also  the  Jacobian  will  illustrate  the  sparse 
structure  shown  in  equation  3.3.  Finally  a  brief  comparison  of  the 
performance  of  algorithms  based  on  each  approach  will  be  presented. 
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Recall,  then,  the  following  CCM  equations  from  the  example 
presented  in  Chapter  2 


In  addition  to  the  information  provided  in  equation  4.1  the  tableau 
approach  requires  a  right  inverse  for  and  a  basis  for  the  null  space 
of  These  can  be  computed  very  efficiently  using  singular  value 
decomposition  techniques  Cl 53  such  as  those  found  n  IMSL  [16]  or 
LINPACKC173.  In  this  example  is  small  enough  to  compute  the 
following  by  inspection. 
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If  the  above  information  is  substituted  into  the  fault  diagram 
equations  (equation  3.1)  the  following  set  of  non-linear  simultaneous 
equations  result: 


45 


1 


JF<x)  = 


0 

0 

0 

0 

0 


2 

I 

0 

0 

0 

3 

7 


0 

a1(1)-a2<1>+  | 

0 

0 

0 


<2)-02(2)+  ^ 
2 


0 

0 


0 

0 

l 

3 

0 

0 

0 

1 

7 

0 


0 

0 

0 

1 

0 

0 

0 

1 


(4.7) 


The  partition  on  the  matrix  in  equation  4.7  emphasizes  the  sparse 
structure  for  the  Jacobian  presented  in  equation  3.3  A  potential 
disadvantage  of  the  sparse  matrix  tableau  formulation  should  be  noted  at 
this  point.  Although  the  Jacobian  in  equation  4.7  is  sparse  the  number 
of  non-zero  entries  is  not  small  enough  to  provide  any  obvious  advantage 
over  the  transfer  function  matrix  approach  in  memory  requirements  to 
store  the  Jacobian.  It  remains  to  be  seen  if  this  short  coming  is 
intrinsic  to  the  sparse  tableau  approach  or  if  it  may  be  circumvented  by 
somehow  exploiting  the  regular  structure  of  the  Jacobian  (equation  3.3). 

Although  the  sparse  tableau  formulation  produces  a  set  of  fault 
diagnosis  equations  for  this  example  (equation  4.5)  which  are  equivalent 
to  the  equations  generated  by  the  composite  transfer  function  matrix 
approach  (equation  5.6  of  Chapter  3)  there  are  several  important 
differences  which  give  the  tableau  approach  numerical  advantages.  First 
the  sparse  tableau  formulation  does  not  require  the  computation  of  a 
transfer  function  matrix  and  hence  the  equations  describing  the  circuit 


are  Much  easier  to  develop.  In  particular  the  inverse  of  the  matrix 
<I-Z(s/r)L^)  ^  is  not  needed.  Second  the  polynomial  order  of  the  fault 
diagnosis  equations  in  the  sparse  tableau  formulations  does  not  increase 
with  the  size  of  the  system  as  the  equations  in  the  composite  transfer 
function  matrix  approach  do.  In  many  applications  (e.g.  RLC  circuits) 
the  sparse  tableau  form  of  the  fault  diagnosis  equations  will  be 
quadratic. 

One  final  comparison  which  should  be  noted  is  the  performance  of 
the  Newton  Raphsons  iteration  scheme  as  applied  to  the  two  different 
approaches.  The  basis  for  the  comparison  is  two  FORTRAN  programs 
executed  on  a  VAX  11/780  minicomputer.  The  first  program  uses  the 
Transfer  Function  matrix  approach  and  the  second  uses  the  Sparse  Tableau 
approach.  Both  approaches  were  applied  to  the  example  presented  in  this 
and  the  previous  chapter  except  that  a  single  complex  test  frequency  (s 
*  jl)  was  used  instead  of  two  real  test  frequencies.  (A  discussion  on 
the  use  of  complex  test  frequencies  will  be  presented  in  the  next 
chapter.)  Table  lisa  summary  of  the  comparison  which  lists  the  number 
of  iterations  required  by  each  scheme  to  find  the  given  parameter  vector 
using  the  nominal  value  (all  parameters  equal  to  unity)  as  a  first  guess 
for  the  Newton  Raphson  iteration. 

The  information  in  Table  1  indicates  that  for  the  example  used  the 
Tableau  equations  converge  faster.  In  on«  case  the  Transfer  Function 
Matrix  equations  failed  to  converge.  These  results  certainly  indicate 
that  further  study  of  the  numerical  properties  of  the  tableau  fault 
diagnosis  equations  is  warranted. 


rcult  to  Illustrate  the  use  of  the  Fault  Diagnosis  Equations. 
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5.  Summary  and  Conclusions 

The  sparse  tableau  formulation  produces  a  set  of  fault  diagnosis 
equations  equivalent  to  those  of  the  composite  transfer  function  matrix 
approach  but  they  are  much  easier  to  compute  and  potentially  easier  to 
solve.  The  development  of  the  Sparse  Tableau  approach  however  is  far 
from  complete  and  before  any  final  conclusions  are  made  on  its  merits 
many  aspects  of  the  approach  remain  to  be  investigated.  Among  these 
are: 

(1)  The  extension  of  the  measure  of  testability  from  the  transfer 
function  matrix  approach  to  the  tableau  approach; 

(2)  The  investigation  of  test  frequency  selection  in  the  sparse 
tableau  context; 

(3)  The  development  of  a  method  for  test  point  selection; 

(4)  A  more  detailed  comparison  of  the  solution  algorithms  of  the 
two  approaches; 

(5)  The  investigation  of  the  modification  of  the  Sparse  Tableau 
Solution  Algorithm  to  exploit  the  quadratic  order  which 
characterizes  the  tableau  equations  in  many  cases; 

(6)  Application  of  the  sparse  tableau  equations  to  the  diagnosis 
of  a  relatively  large  example  (more  than  30  components) 

The  first  item  and  some  aspects  of  the  second  form  a  major  part  of  the 
next  chapter  of  this  report. 
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CHAPTER  4 

CURRENT  RESEARCH  IN  THE  TABLEAU  APPROACH 


1.  Introduction 

This  chapter  contains  a  collection  of  results  stemming  from  recent 
research  into  the  Tableau  Approach  to  fault  diagnosis  [93.  These 
include  a  discussion  of  the  Newton  Raphson  iteration  applied  to 
overspecified  equations/  a  test  for  diagnosabi lity  in  the  Tableau 
context/  bounds  on  test  frequencies  and  the  presentation  of  several 
examples. 


2.  Newton-Raphson  Iterations  and  the  Non-Square  Jacobian 

In  all  previous  discussions  of  the  solution  of  the  nonlinear 
equation  F(x)  =  ®  the  underlying  assumption  was  that  F(x>  mapped  Rn 
into  Rn.  In  this  case  the  Newton  Raphson  iteration  step 


.  ,  <k)w  <k+1  >  <k). 

Jp(X  )<x  -x  )  = 


.  ,  ,  <k+1) 

has  a  unique  solution  for  <x 

non-singular.  In  many  cases  the 

mapping  from  Rn  into  Rm  where  m*n. 


-  F(x<k>)  (2.1) 

x*k*)  when  the  Jacobian  Jp(x*k>)  is 
fault  diagnosis  equation  will  be  a 
To  see  why  this  occurs  recall  the 


tableau  equations  from  Chapter  3: 
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(v)  B.  =  b0(s.) 

In  the  case  of  real  test  frequencies  (complex  test  frequencies  are  dis¬ 
cussed  later  In  the  chapter),  F(x):  DC  RN+pq  *  RMq  and 

jF(Xk>  c  R<"q>x(N*pq> 

where  (1)  N  *  dim  (r) 

<ii)  p  *  dim  (a.-),  i  =  1,2,...,q 
(iii)  M  *  row  dimension  of  Z(s.,r) 

(iv)  q  =  the  number  of  test  frequencies 
Since  N,  p,  M  and  q  are  all  Integers  it  is  not  necessarily  true  that: 

N  +  pq  *  Mq  (2.3) 

If  equation  2.3  holds  the  solution  of  equation  2.1  proceeds  using  any  of 

the  many  accurate  numerical  techniques  available  for  inverting  the 

(k) 

square  non-singular  Jacobian,  Jf(x  ). 
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Under  the  diagnosabi lity  criterion  to  be  developed  later  in  the 

chapter,  a  system  for  which  Jp(»)  does  not  have  full  column  rank  is  not 

guaranteed  to  be  diagnosable  and  therefore  the  usual  techniques  are  not 

directly  applicable.  This  situation  is  characteristic  of  a  lack  of 

sufficient  measurement  data  for  diagnosis.  For  the  underspecified  case 

it  is  necessary  to  augment  the  available  data  by  adding  test 

frequencies,  providing  more  test  outputs,  or  assigning  known  or  assumed 

(k) 

values  to  a  subset  of  the  parameters.  Assuming  however,  that  Jp(x  ) 
has  full  column  rank  the  only  other  possible  relation  for  the  dimensions 
of  JF(x(k>)  is 

N+pq  <  Mq  (2. A) 

A  simple  example  will  illustrate  the  nature  of  the  Newton  Raphson 
iteration  when  inequality  2.4  holds  (Jp(*)  has  more  rows  than  columns). 
Suppose  F(x):  D  C  ♦  R^  is  defined  as 

(2.5) 

F (x)  *  e,  x  =  (x^x^)*. 


F(x^,x^)  = 


*1  +  *2  “  2 


2x1  +  x2  -  3 
x*  ♦  x*  -  2 


The  point  (1,1)*  is  a  solution  to  the  equation 

3x2 

The  Jacobian,  Jp(x^,x2)  e  R  is: 


* 
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JF<X1'X2>  " 


1 


1 

2x. 

i 

2x- 


1 


(2.6) 


Let  the  starting  point  for  the  Newton  Raphson  iteration,  be  (2,2)*. 
Equation  2.1  evaluated  at  the  first  iteration  of  this  example  becomes: 


'4  r 

2  4 

4  4 


X 

5 

6 


(2.7) 


where 


V 

r,«> .  2i 

1  c 

h. 

s 

(i)  , 

2 

x2  -  2 

»  m 

. 

Since  the  vector  -[  4  5  63*  is  not  contained  in  the  space  spanned  by 

the  "column  vectors"  [4  2  43*  and  Cl  4  43*  equation  2.7  is  incon¬ 
sistent. 

How  then  is  the  Newton-Raphson  iteration  step  solved  when  the 
linearized  equations  are  inconsistent?  A  reasonable  approach  to  provide 
a  solution  to  equation  2.1  when  it  is  inconsistent  is  to  find 
x  -  x  =  h  such  that 


||JF(x(k))h  +  F  (xCk>) 1 12 


(2.8) 


is  minimized,  where  | | * | | 2  is  the  Euclidean  norm.  This  "least  squares" 
problem  is  readily  solved  using  singular  value  decomposition  techniques 
C153  available  in  many  computer  subroutine  libraries  such  as 
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LINPACKC17D.  To  see  why  such  a  scheme  is  reasonable  suppose  F(»): 
0  C  Rn  -►  Rm  and  m  <  n.  Consider  the  Taylor  expansion  of  F(x*)  about  the 
point  x  where  F <x*>  =  e. 

F(x*)  =  8  *  F(x)  ♦  JplxXx*  -  x)  ♦  e  (2.9) 

The  term  c  represents  the  effect  all  the  higher  order  terms  of  the 

expansion.  Recall  from  the  derivation  of  the  Newton  Raphson  scheme 

(Section  4  of  Chapter  2)  that  c  is  assumed  to  be  small  enough  to  be 

disregarded.  Before  doing  that  decompose  F(x)  and  e  into  orthogonal 

components  in  the  following  manner: 

F (x)  =  F(x)C  +  F(x)1  (2.10a) 

c  =  eC  ♦  e1  (2.10b) 

where  the  superscript  C  denotes  the  projection  onto  the  column  space  of 
Jp(x)  and  the  superscript  1  denote  the  projection  onto  the  orthogonal 
complement  of  the  column  space  of  Jp(x).  Substituting  equation  2.10 
into  2.9  yields: 

e  =  F(x)C  ♦  F(x)1  ♦  Jp(x*  -  x)  +  eC  +  e1  (2.11) 

Since  those  terms  which  lie  in  the  column  space  of  Jp(x)  are  orthogonal 
to  those  which  lie  in  its  orthogonal  complement,  equation  2.11  can  be 
decomposed  into  two  equations: 

e  a  F (x)^  +  Jp(x*  -  x)  ♦ 


(2.12a) 
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6  =  F(X)1  +  eA 


(2.12b) 


Cl  c  1 

Since  c  =  c  +  c  where  e  and  e  are  orthogonal  then 
Hell*  =  I |eC| ||  +  | |eA| l|.  If  ||e|l2  is  small  then  ||eC||2  and  ||eA|l2 
must  also  be  small.  Neglecting  the  term  eC  in  equation  2.12a  leads  to 
the  following  iteration  step: 


JF(x<k))Cx(k+1)  -x<k>)  = 


-  F(x<k>)C 


(2.13) 


Equation  2.13  is  simply  the  Newton  Raphson  step  used  earlier  except  the 

(k) 

right  hand  consists  of  only  that  part  of  -F(x  )  which  lies  in  the 

(k) 

column  space  of  Jp(x  ).  Not  only  is  equation  2.13  consistent  but  it 
is  mathematically  equivalent  to  the  "least  square"  problem  of  equation 
2.8. 

One  would  expect  that  equation  2.13  has  convergence  properties 

similar  to  the  Newton  Raphson  iteration  when  F:  0  c  Rn  ♦  Rn.  To 

illustrate  this  a  simple  "attraction  theorem”  for  equation  2.13  will  be 

presented  later  in  this  chapter.  Intuitively,  however,  as  the  estimate 
(k) 

x  approaches  the  solution  point  the  linear  term  of  equation  2.9 

becomes  more  accurate  causing  ||e|l2  t0  decrease.  This  implies  that 

i  c  1C 

lie  I l2  and  ||e  Il2  must  also  decrease  since  e  and  c  are  orthogonal 

components  of  e. 


3.  Piagnosabi lity 

The  purpose  of  this  section  is  to  develop  a  test  to  determine  under 
what  conditions  it  is  reasonable  to  attempt  to  solve  the  Tableau  fault 
diagnosis  equations  (Equation  2.2).  Unlike  the  Matrix  Transfer  Function 
Equations  the  Tableau  equations  have  the  ambiguity  unknowns,  a^. 
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131.2.. .,q  in  addition  to  the  parameter  vector,  r.  Because  the 
objective  is  to  identify  the  parameter  vector  r,  computing  then's  may 
seem  unnecessary.  The  role  of  the  these  auxiliary  unknowns  is  to 

U 

"unfold"  the  fault  diagnosis  equations  from  R  where  r  lives  (as  in 
other  approaches)  to  the  larger  dimensional  space,  RN+Pq,  where  x  lives. 
In  RN  the  fault  diagnosis  equations  are  highly  nonlinear.  In  RN+pq  the 
equations  are  often  quadratic,  and  the  data  of  the  end  of  Chapter  3 
supports  the  conjecture  that  solving  for  the  larger  dimensional  x  is 
preferable  to  solving  for  r  directly.  See  Table  1  of  Chapter  3,  page 
49. 

Like  the  fault  diagnosis  equations  in  the  Transfer  Function  Matrix 
Context  (equation  2.7  of  Chapter  2)  the  solvability  of  the  tableau 
equations  (equation  2.2)  is  based  on  the  existence  of  an  isolated 
solution  (Section  3  of  Chapter  2).  For  the  tableau  equations  however 
the  issue  is  confused  by  the  presence  of  the  auxiliary  unknowns,  i  = 

1.2.. ..,q.  For  example  suppose  x  e  R  ,  x  =  (.aj '.Sgr*  r  )  is  a 
solution  to  equation  2.2.  Suppose  furthermore  that  for  every  6  c  R, 
fi  >  0,  there  exists  a  vector  h  e  RN+pq,  with  | | h | | ^  >  0,  such  that 
x*  +  fih  is  also  a  solution.  Clearly  in  this  situation  x*  is  not  an 
isolated  solution.  Now  suppose  the  vector  h  in  the  above  situation 
always  has  the  form  h  =  (h^,h2,...,hq,0).  In  other  words  even  though 
the  solution  x*  is  not  uniquely  specified  by  equation  2.2  the  subvector, 
r,  is  uniquely  specified.  If  such  a  situation  were  possible  a  criterion 
of  solvability,  which  required  that  the  solution  x*  be  isolated  in  order 
to  completely  specify  r,  might  be  too  restrictive.  The  objective  of  the 
next  few  paragraphs  is  to  show  that  such  a  situation  will  not  occur  and 
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thus  permit  the  criterion  of  an  isolated  solution  to  be  used  as  the 

basis  for  a  definition  of  diagnosabi lity. 

Recall  that  the  tableau  fault  diagnosis  equations  are  based  on 
M 

measurements  y  (s.)  in  response  to  inputs  u(s^),  i  *  1,2,...,q.  Using 
the  transfer  function  matrix, 

S(s,r>  =  L22  +  L21 [l-Z<s/r)L11]"1  Z(s,r)L12  (3.1) 

the  measurement  can  be  expressed  in  terms  of  the  parameter  vector 
directly  as: 

y"<s.)  =  [l22  +  L21Cl-Z(sKr)L11  D-1  Z(s,r)L12]  ufs^  (3.2) 

for  i  =  1,2,...,q.  The  family  of  equations  in  equation  3.2  share  the 
same  properties  as  the  fault  diagnosis  equations  of  the  Transfer 
Function  Matrix  approach.  In  fact  the  input  u(s.)  and  s^  i  =  1,2,...,q 
can  be  chosen  so  that  equations  3.2  correspond  exactly  to  the  Transfer 
Function  Matrix  fault  diagnosis  equations  (equation  2.6  of  Chapter  2). 

The  concern  that  basing  a  criterion  for  solvability  of  the  Tableau 
Fault  diagnosis  equation  (equation  2.2)  on  the  existence  of  an  isolated 
solution  can  now  be  dispelled  with  the  following  theorem: 

Theorem  3.3:  The  vector  r*  is  an  isolated  solution  of  equations  3.2  if 

and  only  if  there  exist  auxiliary  vectors  i  *  1,2,...,q  such  that 
*  *  *  *  *  t 

x  =  \®i 'Jli ' •  •  •  * jiq' r  '  1S  *n  isolated  solution  of  equation  2.2.  The 
proof  of  theorem  1  requires  the  following  Lemmas. 

Lemma  3.4:  If  r*  is  a  solution  of  equation  3.2  then  there  exist  unique 
i  *  1^2, ...,q,  such  that  x  =  (^#^^#t.. r  )  is  a  solution  to 
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equation  2.2. 

Proof  of  Lemma  3.4:  The  proof  is  by  construction.  Since  r*  is  given  it 
may  be  used  in  the  tableau  equation  (equation  2.1a  of  Chapter  3)  to 
simulate  the  system.  Now  suppose  that  the  test  frequencies,  s., 
i  *  1,2,...,q,  are  neither  circuit/system  poles  nor  component  poles.  In 
this  case  equation  2.1a  of  Chapter  3  is  well  defined  and  can  be  solved 
uniquely  for  b(s.),  i  =  1,2,...,q.  Denote  these  solutions  b*(s.), 
i  =  1,2,...,q.  From  equations  2.8  of  Chapter  3: 

v£i  =  b*<s.)  -  bg(s.)  i  *  1,2,. ..,q.  (3.5) 

By  construction  the  right  hand  side  lies  only  in  the  column  space  of  V. 
This  implies  that 

RankCVD  =  Rankjv  |  b*^)  -  bQ(s.)J  (3.6) 

and  therefore  solutions  i  *  1,2,...,q,  of  equation  3.5  exist.  Since 
V  has  full  column  rank  these  solutions  are  unique.  Furthermore  the 
construction  insures  that  the  vector  x*  =  satisfies 

equation  2.2. 

Q.E.D. 

Lemma  3.7;  If  x*  =  a  so*-ut*on  t0  equation  2.2  then 

r*  is  a  solution  to  equation  3.2. 

Proof  of  Lemma  :5 .7: 

it  if 

Since  x  is  a  solution  to  equation  2.2,  the  vectors  r  and 
i  *  1,2,...,q,  satisfy 


i 
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[z^r*) 


-»] 


L11Vii  +  W 


s-i 


-  b0(s.)  =  e 


(3.8) 


Recall  that 


W  =  L2?[yKc*l)  -  L22u(i}] 


(3.9) 


and 


W  =  L1lW  +  L12u<si)  (3.10) 

Multiplying  out  equation  3.8  and  substituting  equation  3.10  yields: 

Z(s.,r*)L11[b0(si)  ♦  Vo*]  +  Z(s.,r*)L12u(s.)  -  |b0(s.)  ♦  «*J  =  e(3.11) 
Solve  for  bQ(s.j)  +  Vou : 

b0(si)  +  v“*  =  [l-Z(si,r*)L11]'1Z(s.,r*)L12u(si)  (3.12) 

Multiply  on  the  left  by  L?1  and  recall  that  the  columns  of  V  span  the 
null  space  of  L21.  Hence 

L2lW  =  L2i[I"z(si'r*>l-11]"1z<sirr*)L12u(s.)  (3.13) 

Finally  substitute  equation  3.9  into  3.13  to  obtain: 

y"(s.)  =  <L22  +  L21[l-Z(si,r*)L11J'1Z(si<,r*)L12>  u(s.)  (3.14) 

i  *  1,2,...,q  for  which  coincides  with  equation  3.2. 


Q.E.D 
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Proof  of  Theorem  3^:$: 

Part  1_  "  "only 

*  *  *  ★  * 

Suppose  x  =  'Sb\'!i2'mmm'SLq'r  ^  1S  n°*  an  isolated  solution  of 
equation  2.2  but  r*  an  isolated  solution  of  equation  3.2.  Since  x* 
is  not  isolated,  for  any  arbitrarily  small  6  >  0  there  exists  another 
solution  8  =  (aJj,a2,...,5q,f>t  such  that  0  <  1 1 x*— S 1 1 2  <  fi'  Now 
according  to  Lemma  3.7  7  must  also  be  a  solution  to  equation  3.2  and 
furthermore  r*  4  P  (If  r*  =  P  then  by  Lemma  3.4,  j**  =  £?,  i  =  1,2,...,q 
which  contradicts  x*  *  8.).  Since  r*  and  P  are  corresponding  components 
of  x*  and  8  respectively  then  0  <  | | x*— 8 1 J ^  <  6  anc'  r*  *  P  imply  that 
0  <  I  |  r*  -  Ml-  <  6*  But  6  can  be  chosen  arbitrarity  small  which 
implies  P  becomes  arbitrarily  close  to  r*  implying  r*  is  not  an  isolated 
solution  to  equation  3.2.  This  contradicts  the  assumption  that  r*  is 
isolated. 

Part  2  -  "vf" 

Now  suppose  r*  is  not  an  isolated  solution  of  equation  3.2.  Use  r* 

in  Lemma  3.4  to  construct  a  solution  to  equation  2.2.  Denote  this 

...  *  .  *  *  *  *.  t  .  .  *.  . 

solution  x  =  (a^/Oj^'-wa^r  )  and  assume  that  x  21  an  isolated 

solution  to  equation  2.2.  From  the  construction  in  Lemma  3.4  the 

i  =  1,2,...,q  and  hence  x  =  (a.^.S2*-«-/<iq,r>  are  continuous  functions  of 

r.  Now  select  an  P  and  construct  the  corresponding  8.  From  the 

property  of  continuity  for  any  e  >  0  there  exists  a  6  >  0  such  that 

whenever  ||r*-P|)2  <  «,  ||x*-S|l2  <  e.  Since  r*  is  not  an  isolated 

solution  the  P  can  always  be  chosen  to  be  a  solution  as  well. 

Consequently  from  Lemma  3.4  8  must  be  a  solution  to  equation  2.2.  Since 

c  is  arbitrary,  in  any  open  neighborhood  of  x*  it  is  possible  to  find 
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another  solution  8.  This  contradicts  the  assumption  that  x*  is 
isolated. 

Q.E.D. 

Lemma  3.4  and  3.7  and  Theorem  3.3  clearly  establish  the  equivalence 
of  the  Fault  Diagnosis  equations  in  the  Tableau  Form  with  those  in 
transfer  function  matrix  form  and  make  the  following  definition 
possible. 

Definition  3.15:  A  circuit/system  modeled  by  equation  2.1  of  Chapter  3 
is  said  to  be  "DIAGNOSABLE"  if  and  only  if  there  exist  real  test 
frequencies/  s^,  and  corresponding  test  inputs  u(s.j),  i  =  1,2,...,q, 
with  q  finite,  for  which  a  solution  to  equation  2.2  exists  and  is  an 
isolated  solution. 

With  respect  to  the  above  definition  the  terms  DIAGNOSABLE  and 
DIAGNOS ABILITY  will  be  used  in  the  context  of  the  tableau  approach  to 
avoid  confusion  with  the  measure  of  testability  defined  in  the  transfer 
function  matrix  context.  Also  the  test  frequencies,  s^,  are  not 
necessarily  distinct.  For  example  a  system  with  two  inputs  could  be 
diagnosable  with  q=2  where  s^  *  5,  u(s^)  =  (1,0)*  and  s2  *  5, 
u<s2)  =  (0,1)*. 

The  next  theorem  provides  an  explicit  test  for  diagnosabi lity. 

Theorem  3.16:  Suppose  a  circuit/system  is  characterized  by  a  parameter 
vector  r.  Let  x  be  a  solution  to  equation  2.2  which  includes  the 
subvector,  r.  Then  the  circuit/system  is  diagnosable  for  the  parameter 
vector,  r,  if  there  are  sufficient  test  frequencies,  s.,  and 
corresponding  inputs,  uls^),  i  =  1,2,. ..,q,  such  that  Jp(x)  has  full 


column  rank 
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Proof:  Recall  that  in  the  tableau  fault  diagnosis  equation,  F(x)  =  0, 

the  nonlinear  function  F:  0  C  RN+pq  +  RMq  where  Mq  >  N+pq.  Since  the 

Mqx(N+pq)  Jacobian  matrix  has  full  column  rank  it  must  have  an 

(N+pq)x(N+pq>  submatrix,  J  x(s>,  which  is  non-singular.  Define  F: 

D  C  RN+pq  ♦  R^+pq  to  be  a  subvector  of  F  whose  entries  are  the  entries 

of  F  corresponding  to  the  rows  of  J  (x).  Clearly  if  x  satisfies 

F 

F(x)  =  0  it  must  satisfy  F(x>  *  0.  Since  J  (s)  is  non-singular  the  In- 

IE 

verse  Function  Theorem  insures  that  x  must  be  an  isolated  solution  of 
F(x)  =  0.  The  vector  x  must  likewise  be  an  isolated  solution  of 
FCx)  *  0  and  therefore  definition  3.15  is  satisfied.  v 

Q.E.D. 

To  be  of  any  practical  use  the  fact  that  J p ( • >  has  full  column  rank 
at  some  specific  solution  point  must  infer  that  it  has  full  column  rank 
elsewhere.  This  is  important  for  two  reasons.  First,  the  Newton- 
Raphson  iteration  scheme  requires  that  Jp(x)  have  full  column  rank  to 
uniquely  solve  the  Newton-Raphson  step  for  a  sequence  of  estimate  x  . 
If  at  any  point  there  is  no  unique  solution  the  iteration  scheme  could 
not  proceed.  Second,  it  would  be  reassuring  to  know  that  when  a 
circuit/system  is  designed  to  be  diagnosable  for  some  nominal  parameter 
vector  r,  it  will  still  be  diagnosable  should  the  parameter  vector 
change  due  to  aging  or  failure  of  the  components. 

In  practice  the  use  of  the  Newton-Raphson  iteration  scheme  requires 
that  JF(x)  have  full  column  rank  almost  everywhere.  In  this  context  al¬ 
most  everywhere  means  everywhere  except  in  a  lower  dimensional  algebraic 
variety  C93.  If  this  is  true  then  whenever  the  algorithm  encounters  a 
"bad"  value  for  x  an  arbitrarily  small  perturbation  of  x 
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will  yield  a  value  from  which  the  iteration  can  proceed.  The  following 
theorem  shows  that  this  is  indeed  the  case: 

Theorem  3.17:  Suppose  x*  e  RN+pq  and  Jp(x*)  has  full  column  rank  then 
JF(x)  has  full  column  rank  for  almost  all  x  e  RN+pq. 

Proof:  If  JF<x*>  has  full  column  rank  then  there  is  an  N+pq  by  N+pq 
submatrix,  M(x*),  such  that  detCM(x*)D  *  0.  Since  every  entry  of  Jp(x*) 
is  either  a  constant  or  a  linear  combination  of  the  elements  of  x*,  this 
determinant  is  a  polynomial  in  x*.  Now  consider  the  polynomial 
detCH(x)].  It  is  not  zero  at  x  =  x*  and  therefore  is  not  identically 
zero.  The  zero  set  of  this  polynomial  is  an  algebraic  variety  of  RN+pq 
of  dimension  no  larger  than  N+pq-1  and  contains  all  points  where  Jp(x) 
is  NOT  full  column  rank. 

Q.E.D. 

The  question  still  remains:  Can  the  parameter  vector  change  to  make 
a  diagnosable  circuit  non-diagnosable?  Diagnosabi lity  is  defined  only 
for  those  points  x  e  RN+pq  which  are  possible  solutions  to  equation  2.2 
and  according  to  lemmas  3.4  and  3.7  this  is  an  N  dimensional  subspace  of 

RN+pq.  c(.ear(.y  measure  over  RN,pq  is  not  appropriate  to  this  situation. 

N  M 

Since  the  solution  space  is  isomorphic  to  R  the  measure  over  R  is  used 
in  the  following: 

Theorem  3.18:  Suppose  a  circuit/system  is  known  to  be  diagnosable  via 
Theorem  3.16  for  parameter  vector  r*.  Then  it  is  diagnosable  for  almost 
all  r  e  R  . 

Proof:  By  Theorem  3.16  test  frequencies  exist  such  that  Jp(x*)  has  full 
column  rank.  (Note:  x*  is  determined  uniquely  from  r*  as  per  Lemma 
3.4.)  If  Jp(x)  is  evaluated  only  over  the  solution  space  of  equation  2.2 


it  may  be  considered  a  function  of  r  only,  since  the  _a.,  i  s 
are  functions  of  r.  Jp  can  be  expressed  explicitly  as  a  function  of  r 
by  replacing  the  term  CL^V^.  +  ag(s^)D  in  expression  3.4b  of  Chapter  3 
as  follows: 

IV*  *  w]  *  ['  -  Hiz(Vr>]  Li2"<V  <3-1,> 


Now  all  entries  of  Jp(r)  are  constants,  polynomials  in  r,  or  ratios  of 
polynomials  in  r.  Because  Jp(r*>  has  full  column  rank,  k,  there  is  a  k 
by  k  minor  M(r*>  such  that  det  M(r*>  is  defined  and  is  not  zero.  Now 
consider  the  parameter  vector  r  #  r*.  The  determinant  of  MCr)  can  be 
written  as  a  ratio  of  polynomials  in  r: 


P1<r) 

DetCM(r)D  =  ^ 


(3.20) 


it 

Since  M(r)  is  not  zero  and  defined  at  r  polynomials  p^(»)  and  P2(»>  are 
not  identically  zero.  Jp(r)  has  full  column  rank  everywhere  except 
possibly  over  the  zero  sets  of  p^(»)  or  P2(»>,  which  are  algebraic 
varieties  in  R  and  hence  so  is  their  union. 

Q.E.D. 

The  discussion  thus  far  has  required  real  test  frequencies.  The 
most  likely  test  signal  in  fault  diagnosis  testing,  however,  is  a 
sinusoidal  signal  which  is  equivalent  to  using  a  complex  pair  of  test 
frequencies.  It  is  appropriate  therefore  to  illustrate  how  the  fault 
diagnosis  equations  are  affected  by  the  use  of  complex  test  frequencies. 
When  the  s^  =  ju.,  i  =  1,2,...,q  are  complex,  the  fault  diagnosis 
equations  simply  become  complex  vector  equations  but  since  r  is 


l 


■^1 
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generally  real  It  is  convenient  to  use  real  arithmetic.  To  do  this  it 
is  necessary  to  convert  the  fault  diagnosis  equation  to  an  equivalent 
set  of  real  equations.  Equation  2.2  becomes: 

ReCF(xRft)3  =  e 
ImCF(xRj,)3  =  e 

where  ReC  3  and  ImC  3  denote  real  and  imaginary  parts  respectively.  The 

DVl+2pq  . 


'Re 


'Re 


ReCa3 

ImCa] 

r 


is  defined: 


*1 

■»2 


where  a  = 


(3.21) 


The  Jacobian  Jj.(xRe)  has  the  form: 


JF<xRe> 


aReCF3  aReCF3  aReCFTj 
aReC  •  almTaT  ar 

almCFD  almCFD  aImCF3 
■Z“Jr“ 


(3.22) 


Like  the  Jacobian  in  the  real  test  frequency  case,  JF<*Re>  is  sparse  and 

each  block  has  a  diagonal  or  block  diagonal  structure.  In  light  of  the 

above  discussion  it  is  possible  to  state  the  following  test: 

Corollary  3.23:  (to  Theorem  3.16):  A  circuit  characterized  by  parameter 

vector,  r,  is  diagnosable  if  there  are  complex  test  frequencies, 

and  corresponding  inputs,  u(s^),  i  =  1,2,...,q,  such  that 

J,(xD.J  has  full  column  rank. 

F  Re 


J 
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Q.E.D. 

Proof;  Special  *e  of  Theorem  3.16. 

This  last  test  Is  significant  because  a  single  sinusoidal  input 
embodies  two  test  frequencies,  ±jw,  and  thus  fewer  input-output 
measurements  are  required. 


4.  Test  Frequencies 

The  use  of  Theorem  3.16  or  corollary  3.23  requires  the  selection  of 
appropriate  test  frequencies.  Presently  no  method  is  known  to  make  this 
selection  in  an  optimal  manner.  Nevertheless  it  is  possible  to  place 
some  simple  bounds  of  the  number  of  test  input/frequency  combinations 
necessary  to  satisfy  the  conditions  of  Theorem  3.16  and  its  corollary. 
The  bounds  developed  here  are  applicable  to  complex  test  frequencies  and 
greatly  simplify  the  use  of  corollary  3.23  in  determining 
diagnosabi lity. 

Theorem  4.1  (lower  bound):  Suppose  a  circuit/system  has  the  following 
characteristics: 

i)  N  real  parameters 

ii)  DimCnullCL^DJ  =  p, 

iii)  The  test  frequencies  are  s.  and  test  inputs  are  u(s.),  and 
i  =  1 ,2, . • «,q. 

iv)  The  row  dimension  of  Z(»,*)  is  M 
The  minimum  number  of  test  inputs/frequencies,  q,  required  to  make  the 
circuit/system  diagnosable  is 


8 


N 


2(M-p) 


(4.2) 


where  the  square  bracket  indicates  "least  upper- bounding  integer." 
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Proof;  For  complex  frequencies  it  is  convenient  to  consider  a  set  of 
equivalent  real  scalar  equations  since  the  parameter  vector  r  is  real. 
Thus  q  complex  test  frequencies  result  in  2Mq  equations  in  N+2pq 
unknowns  (the  j*.  may  be  complex  but  r  is  not).  In  the  best  case  all 
equations  are  independent.  Hence  2Mq  N+2pq  which  implies 

8- 

Q.E.D. 

Note:  Implicit  in  the  above  proof  is  the  fact  that  N  must  be 
greater  than  p  in  order  for  q  to  be  positive  and  finite.  Intuitively, 
if  the  dimension  of  ambiguity  at  any  single  test  frequency  is  not 
strictly  less  than  the  number  of  scalar  equations  generated  by  each  test 
frequency,  no  choice  of  test  frequencies  could  possibly  resolve  it. 
Theorem  4.3  (upper  bound):  Consider  the  same  circuit/system  as  Theorem 
4.1.  If  it  is  diagnosable  the  maximum  number  of  test  inputs/frequencies 
necessary  is  N. 

Proof:  Each  uf  the  terms  f . (r)g.. (ok)  -  B.  =  0  in  equation  2.2  consists 
of  N  complex  scalar  equations  in  N+p  unknowns.  N  unknowns  corresponding 
to  the  parameter  subvector  r  are  common  to  each  set.  Converting  to  an 
equivalent  set  of  real  equations  this  yields  2K  equations  in  N+2p 
unknowns  (again  r  is  real).  Therefore  the  first  test  input /frequency 
produces  N+2p  unknowns  and  each  subsequent  measurement  introduces  2p 
more.  Clearly  if  each  such  set  of  2M  equations  contributes  2p 
independent  equations  to  the  entire  set  there  could  be  no  locally  unique 
solution,  whereas  if  each  contributes  2p+1,  a  locally  unique  solution  is 
possible.  Since  the  circuit/system  is  diagnosable  by  hypothesis  then 


the  worst  case  (requiring  the  most  test  frequencies)  is  2p+1  independent 


equations  per  test  frequency.  In  this  case  q  must  satisfy: 
q(2p+1)  £  N+2pq  and  q  £  N.  Since  q=N  is  the  minimum  integer  which 
satisfies  the  worst  case  this  must  be  the  most  test  frequencies 
necessary. 

Q.E.D. 

For  a  given  circuit  is  now  possible  to  choose  test  frequencies  and 
determine  if  the  choice  renders  the  circuit  diagnosable  for  some  nominal 
set  of  parameters.  After  a  circuit  is  found  to  be  diagnosable  the 
Newton-Raphson  iteration  scheme  can  be  used  to  solve  for  the  actual 
values  for  those  parameters  from  measurements  taken  at  the  selected  test 
frequencies. 

5.  Nonlinearity  of  the  Fault  Diagnosis  Equations 

An  extremely  important  special  case  form  of  the  component  transfer 
function  matrix  arises  when 

Z.(s,r.>  =  r.  Z.(s)  (5.1) 

In  this  case  each  of  the  product  terms  of  the  Tableau  form  of  the  fault 
diagnosis  equations  is  affine  C21 3  and  the  components  of  F(x)  are  always 
quadratic  polynomials  in  x,  regardless  of  the  number  of  components  in 
the  system.  This  is  not  the  case  with  the  Transfer  Function  Matrix  form 
of  the  Fault  diagnosis  equations  (equation  2.6  of  Chapter  2).  Each 
entry  of  equation  2.6  of  Chapter  2  has  the  form: 
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* 


k 


Jkn 

pkd 


C5.i 


where  m^  is  the  measurement  of  one  entry  of  the  transfer  function  matrix 
and  P^r)  an<*  pkd^r^  are  P°lyno*'**l*  in  r.  The  ratio  ^(rJ/py^f)  is 
a  scalar  transfer  function  from  the  system  input  to  the  system  output 
associated  with  the  entry  of  the  transfer  function  matrix  to  which  the 
measurement,  m^,  corresponds. 

As  the  size  of  the  system/circuit  increases  the  degree  of  the 
polynomials,  pkn(r)  and  will  generally  increase.  A  simple 
approach  to  see  why  this  occurs  is  to  form  the  directed  graph  which  is 
equivalent  to  the  CCM  equations  and  then  use  Mason's  gain  formula  C20D 
to  make  some  general  observations  about  the  nature  of  p^n<r>  and  p^^< r> 
as  the  number  of  components  of  the  circuit  in  question  increases. 

A  directed  graph  is  easily  constructed  from  the  CCM  equations  by 
assigning  to  each  component  a  directed  branch  with  weight  r-Z^(s), 
i  =  1,2,...,N.  Thus  each  branch  weight  is  a  first  degree  polynomial  in 
r.  The  vertices  of  each  of  these  branches  represent  the  component’s 
input  and  output.  The  graph  is  completed  by  connecting  the  component 
inputs  and  outputs  and  vertices  representing  the  system  inputs  and 
outputs  with  directed  branches  whose  weights  correspond  to  the  non  zero 
entries  of  the  (connection  matrices)  in  such  a  way  that  the 
appropriate  conservation  laws  are  satisfied.  Figure  8  is  an 
illustration  of  a  graph  constructed  in  such  a  manner.  The  graph  is 
based  on  the  CCM  equations  for  the  example  circuit  from  Chapter  2 
(equations  5.1,  5.2  and  5.3  of  Chapter  2). 
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Before  starting  Mason's  gain  formula  the  following  definitions 
associated  with  the  use  of  graphs  are  necessary  C203: 

Path  gain  -  The  product  of  the  weights  of  the  branches  of  a 

directed  path  In  the  graph. 

Cycle  gain  -  The  product  of  the  weights  of  all  branches  in  a 

directed  circuit  of  the  graph. 

Let 

tj  *  sum  of  cycle  gains  of  all  directed  circuits  in  the  graph; 

t^  =  sum  of  products  of  all  vertex-disjoint  (non-touching)  directed 
circuits  tak.-n  two  at  a  time; 

tk  =  sum  of  products  of  all  vertex-disjoint  directed  circuits  taken 
k  at  a  time; 

A  *  1-t1  ♦  t2  -  t3  ♦  •••  ♦  <-1>Q  t 

where  q  is  the  maximum  number  of  vertex-disjoint  directed 
circuits  in  the  graph; 

*  path  gain  of  the  kth  directed  path  from  output  to  input  vertex 
for  which  the  transfer  function  is  to  be  computed; 

=  the  value  of  A  for  that  part  of  the  directed  graph  having  no 
vertices  in  common  with  the  path  associated  with  P^. 

Notice  that  the  expressions  for  A,  A^  and  P^  are  all  polynomials  in  r. 
The  scalar  transfer  function  can  now  be  computed  using  Mason's  gain 
formula  C20D  as: 

pkn<r>  . 
pJjTrJ  *  t 


(5.3) 
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Consider  each  of  the  terns  A,  A^  and  as  the  circuit/system  and 
hence  its  graph  increase  in  size.  The  paths  from  input  vertices  to 
output  vertices  will  tend  to  get  longer  and  the  number  and  size  of  the 
directed  circuits  will  likewise  increase.  Thus  the  polynomial  degree  of 
A,  A^  and  and  consequently  of  p^fr)  and  p^(r)  will  increase  in 
proportion  to  the  size  of  the  circuit/system. 

The  degree  of  the  components  of  the  Tableau  Equation  on  the  other 
hand  remains  fixed  (quadratic)  as  the  system/circuit  size  increases. 
Therefore  each  term  has  a  constant  second  derivative  with  respect  to  r 
and  higher  derivatives  vanish.  A  solution  algorithm  which  exploits  this 
easily  computed  second  derivative  information  should  offer  greatly 
improved  performance  over  the  Newton-Raphson  scheme  discussed  earlier. 


6.  Convergence  of  the  Newton  Raphson  Iteration 

At  the  beginning  of  the  Chapter  a  form  of  the  Newton-Raphson 

iteration  step  (equation  2.13)  was  derived  to  solve  F(x)  =  0,  where  the 

nonlinear  function  F:  0  c  Rn  *  Rm  and  m  >  n.  The  objective  of  this 

section  is  to  show  that  the  iteration  of  equation  2.13  is  locally 

convergent.  An  iteration  is  said  to  be  locally  convergent  if  the 

iteration  will  converge  to  the  solution  point,  x*,  whenever  the  initial 

estimate,  x*°*,  is  sufficiently  close  to  x*.  The  following  definition 

and  theorem  from  C193  are  the  basis  for  achieving  this  objective. 

Definition  Cl 93 i  Let  6:  D  c  Rn  ♦  Rn.  Then  x*  is  a  point  of  attraction 

(lc+1)  (k) 

of  the  iteration,  x  =  G(x  ),  k  =  0,1,...,  if  there  is  an  open 

neighborhood  S  of  x*  such  that  Sc  D  and,  for  any  xQ  e  S,  the  iterates 

(k)  * 

{x  >  all  lie  in  D  and  converge  to  x  . 
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Ostrowski  Theorem  [193:  Suppose  that  G:  D  c  Rn  ♦  Rn  has  a  fixed  point 
x*  e  int(D)  and  is  Frechet  differentiable  at  x*.  If  p(6  (x*>)  <  1  then 
x*  is  a  point  of  attraction  of  the  iteration  x*k+^*  =  G(x^*).  Where 

(i)  x*  is  a  fixed  point  of  G  if  x*  s  G(x*> 

I 

(ii)  G  is  the  Frechet  derivative  of  G. 

(iii)  p(A)  is  the  magnitude  of  the  largest  eigenvalue  of  A  (known  as 
the  spectral  radius). 

To  establish  the  local  convergence  of  equation  2.13,  use  the 
following: 

Theorem  6.1 :  Suppose  there  is  a  point  x*  e  Rn,  such  that  F(x*)  =  8, 
where  F:  0  c  Rn  ♦  Rm  and  m  >  n.  Suppose  also  that  in  a  neighborhood  of 
x*  the  function,  F,  has  a  continuous  Frechet  derivative,  Jf,  and  has 
full  column  rank.  Then  x*  is  a  point  of  attraction  for  the  iteration 

„<K*1>  .  „<k>  .  jJ<x<k)>FCx<k>>  (6.2) 

Proof:  Define  G:  D  c  Rn  ♦  Rn  as 

G(x)  *  x  -  |jJ(x)Jp(x)J  Jp(x)F(x)  (6.3) 

Notice  that  x*  is  a  fixed  point  of  G  -  i.e. 


x*  =  G(x*)  (6.4) 

For  convenience  note  that  a  left  inverse  of  Jf(x)  denoted  CJp(x>]  L  C14D 
has  the  form. 


[jp(«)]  L  *  [jJ(x)JF(x)] 


jJ(x) 


(6.5) 


and  notice  that 
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JF(x)  =  I 


(6.6) 


To  apply  the  Ostrowski  Theorem  it  is  sufficient  to  prove  that  G*(x*)  *  e 
in  which  case  the  spectral  radius  is  zero.  Equivalently  show  that 

lim  - — -  |  |G(x>  -  G(x*)||-  =0  (6.7) 

Ilx-xll2  2 

To  prove  equation  6.7  holds,  evaluate  the  argument  of  the  limit  opera¬ 
tion  using  equation  6.3  and  note  that  F(x*)  =  6.  In  particular 

- ^ -  |  |G(x)-G(x*)  |  U  =  - 1- -  1 1 1  (x-x*)  -  L(x)l  F(x)||,  (6.8) 

l|x-x*|l2  2  ||x-x*|l2  LF  -I  2 

Next  use  equation  6.6  and  the  fact  that,  F(x)  =  F(x)  -  F(x*)  which  pro¬ 
duces 


- - -  |  |G(x)  -  G(x*)  1 1?  * 

1 1 x— x* 1 1 2  2 


x-x 


[JF(X)]  L  [jF(x)(x-x*)  -  [f (x)  -  F(x*)]J 


(6.9) 


and 


- S -  |  I G (x)  -G(x*)||,  < 

1 1  x— x*  1 1 2  2 


Jp(x)(x-x*) 


[F(x)  -  F(x^||2 


(6.10) 
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r  ! 


t  ; 

* 


where 


-L 


max  j|JF(x)“Lh||2 


(6.11) 


is  the  spectral  norm  of  the  matrix  Jf(x)"L.  Add  and  subtract  Jp(x)(x- 

x*)  from  the  last  term  of  equation  6.10  and  use  the  triangle  inequality 
to  obtain: 

1  'Ig(x)  -  G(x*)||-  < 


x-x 


(x)  -  JfCx*)) (x-x*)  2  + 


(6.12) 


I  lx 


Since  Jp(x)  is  continuous 


lim 


x-x*  .  e  1 1*'*  "2 


and  by  definition  of  the  Frechet  derivative 

1 


jL -  | jtJpCx)  -  JF(x*)](x-x*>  2  =  0  (6.13) 


lim 

*  .  1 1 x-x"] |- 

x-x  *  6  1 '  '2 


jp(x*)(x-x*)  -  F(x) 


F(x  ) |  j  s  0  <6.14) 


'  * 

Hence  G  (x  )  =  0. 


( 

J 


Q.E.D. 
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This  establishes  the  validity  of  the  least  squares  solution  for  the 
Newton  Raphson  iteration. 

This  section  concludes  the  theoretical  development  of  the  tableau  ap¬ 
proach  to  date.  The  next  section  will  demonstrate  the  application  of  this 
theory  to  several  examples. 

7.  Fault  Diagnosis  Examples 

The  purpose  of  this  section  is  to  illustrate  the  use  of  the  fault  diag¬ 
nosis  equations  by  applying  them  to  several  example  problems.  Each  example 
includes  its  CCM  characterization  as  well  as  nominal  and  actual  parameter 
values.  Currently  no  method  exists  in  the  context  of  the  tableau  approach 
for  optimal  test  point  selection.  In  each  example  the  selection  of  test 
points  took  place  via  the  trial  and  error  procedure  below: 

(1)  Model  the  ci rcuit/system  in  terms  of  the  CCM. 

Note:  This  includes  the  use  of  frequency  and  impedance  scaling.  As  in  the 
simulation  process  the  scaling  enhances  the  numerical  variations  in  the  nom¬ 
inal  values  of  the  parameters  from  component  to  component.  Such  scaling  is 
always  possible  and  requires  only  that  actual  measurements  undergo  a  reverse 
scaling  before  their  entry  into  the  solution  program. 

(2)  Select  candidate  test  outputs  (usually  functional  outputs). 

(3)  Select  a  minimum  set  of  test  input /frequency  combinations  (via 
Theorem  4.1 ) . 

Note:  All  the  examples  presented  here  are  single  input  circuits/systems. 
Selecting  test  input /frequency  combinations  amounts  to  selecting  distinct 
test  frequencies.  In  all  cases  u(s . )=1i-j0,  i=1,2,...,q. 

(4)  Test  the  rank  of  Jp(x°>  (see  Theorem  3.16)  to  determine  if  the 
circuit/system  is  diagnosable,  where  x°  is  the  nominal  parameter 


vector 
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(5)  If  the  diagnosabi lity  test  fails,  include  an  additional  test 
input /frequency  combination,  if  the  test  succeeds  go  to  step  8. 

(6)  Repeat  steps  4  and  5  until  N  input /frequency  combinations  are  used 
(see  Theorem  4.3). 

(7)  Select  additional  test  outputs  and  return  to  step  3. 

(8)  Use  the  resulting  test  outputs  and  test  input/frequency  combinations 
to  generate  the  tableau  fault  diagnosis  equations  (equation  2.2). 

The  solution  of  the  resulting  equations  proceeds  via  a  computer  imple¬ 
mentation  of  the  Newton-Raphson  iteration  step  with  the  nominal  parameter 
values  and  corresponding  ambiguity  vectors  (a^,  i=1,...,q)  serving  as  the 
first  estimate  of  the  solution.  This  "diagnosis  program",  written  in  FOR¬ 
TRAN  uses  the  subroutine  LLSQF  from  INSLC163  to  determine  the  least  squares 
solution  to  the  Newton-Raphson  iteration  (equation  2.1).  Sparse  matrix  tech¬ 
niques  were  not  utilized.  The  program  executed  on  a  VAX  11/780  computer. 
Each  example  concludes  with  a  summary  of  the  performance  of  this  solution 
algorithm.  The  next  chapter  introduces  a  modification  to  the  solution  algo¬ 
rithm  which  exploits  the  quadratic  form  of  the  tableau  fault  diagnosis  equa¬ 
tions  and  several  of  the  examples  presented  here  will  also  demonstrate  the 
performance  of  the  modified  algorithm. 

Example  7.1;  The  first  example  is  based  on  the  AC  equivalent  circuit  of 
the  single  stage  transistor  amplifier  in  figure  9  selected  from  reference 
C63.  The  component  connection  equations  are  shown  below:  (Note:  The  input 
voltage  for  the  dependent  source  is  defined  differently  thsn  it  was  in  C63 
making  Z(s,r)  diagonal  and  changing  one  entry  in  L^.) 


Figure  9.  Amplifier  model  for  example  7.1. 
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A  basis  for  the  null  space  of  L-^  in  matrix  form  is: 


For  this  circuit  M=N=12  ,  p=9  and  q®2.  ALL  parameters  have  nominal  values 
of  unity.  For  test  frequencies  s1  ■  jl .4  and  s2  *  j2  the  nominal  values  for 
the  ambiguity  vectors,  and  are: 
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■  .72545 E-01  ‘ 

-.42669E+0(f 

.26209E+00 

.28641  E-01 

.92251 E-01 

.28133E-01 

.18192E+00 

-.38153E-01 

“1 " 

-.39387 E-01 

♦  ) 

.12915E+0Q 

.64166E-01 

-.98703E-01 

.92251 E-01 

.281 33 E-01 

•16710E+00 

.11155E+00 

.261 66 E+00 

.  .21 033  E+00. 

-.86316E-01  ' 

-.48784E+0(f 

. 29425 E+00 

.80346 E-01 

.10092E+00 

.10641 E-01 

•17635E+00 

-.11167E-01 

-.21 282 E-01 

♦  j 

.201 84 E+00 

.35081 E-01 

-.13916E+00 

.10092E+00 

.10641 E-01 

.19156E+00 

.97295E-01 

.31 972  E+00 

.21512E+00 

Next  the  parameter  vector  r  was  changed  to: 

r  *  (2.  1.2  .9  .9  1.1  .8  1.  1.1  2.  .9  1.2  1.1)t  (7.7) 


The  small  variations  from  unity  represent  normal  production  variations  for 
"good"  components  while  the  larger  variations  for  the  first  and  ninth  com¬ 
ponents  represent  failures.  This  perturbed  parameter  vector  r  (representing 
an  actual  circuit  under  test)  was  used  to  compute  the  following  test  out¬ 
puts: 


-.18299E-01" 

'  . 44648 E-01" 

.24554E+00 

+  j 

•55346E+00 

_  .73848E-01. 

.15117E+00 

(7.8a) 
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y"(s2) 


.6214ZE-0?]  f  .  67445 E-Of 

.41992E+00  ♦  )  . 68044 E+OO 

.10101  E+OOj  L  .19389E+00. 


(7.8b) 


These  test  outputs  and  the  nominal  data  were  the  input  to  the  diagnosis  pro- 
gran  for  determination  of  the  actual  parameter  values.  The  program  used 
twelve  Newton-Raphson  iterations  and  31.5  seconds  on  the  VAX  11/780  to  com¬ 
pute  the  estimates  of  the  circuit  parameters,  all  within  .11%  of  the  actual 
values. 

Exam£l£  7^.£:  This  example  consists  of  the  op-amp  circuit  shown  in  Fig¬ 
ure  10.  The  op-amp/feedback  resistor  combination  forms  a  "voltage-follower" 
which  is  extremely  insensitive  to  changes  in  the  actual  value  of  the  feed¬ 
back  resistor.  Intuitively,  the  individual  diagnosis  of  this  resistor's 
value  will  be  numerically  unstable.  It  is  therefore  reasonable  to  combine 
the  op-amp  and  feedback  resistor  as  a  single  component  for  the  purposes  of 
fault  diagnosis  and  model  them  as  a  voltage  controlled  voltage  source.  The 
resulting  circuit  with  parameters  scaled,  appears  in  Figure  11.  The  test 
points  for  this  example  are  a  result  of  the  trial  and  error  process  dis- 

•D 

cussed  at  the  beginning  of  the  section.  The  CCM  equations  as  well  as  l_21 
and  V  are  given  below: 


■  - 

“  - 

•  - 

b1 

r1 

a1 

b2 

r2s  0 

*2 

b3 

S 

r3 

*3 

b4 

0  r4/s 

a4 

b5 

r5 

a5 

(7.10) 
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0 

-1 

-1 

0 

0 

0 

0 


1 

0 

0 

0 

0 

0 


1 

0 

0 

1 

0 

0 


1  1 


0 

0 

-1 

0 

1 

0 

0 


-1 


n 

V 

1 

b2 

1 

b3 

+ 

0 

b4 

0 

0 

J 

b5 

0 

(7.11) 


L_R 

L21 


m 

0 

0 

0 

1 


0" 

.5 

.5 

0 

0 


(7.12) 


V  = 


"1  0  (51 

0-10 
0-10 
0  0  1 

0  0  0 


(7.13) 


For  this  example  M=N=5,  p=3  and  q=2.  The  test  frequencies  are 

s^  =  j  .5  and  «  j  1.5.  As  in  the  previous  example  the  parameter  vector, 

r*,  used  to  compute  test  outputs  was  a  perturbation  of  the  nominal,  rQ.  The 

values  used  for  r*  and  r  are 

o 


rQ  *  (1  1.5  .5  2.5  1)' 


(7.14) 
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r*  =  <2.2  .75  .90  2.2  I)1 


(7.15) 


The  nominal  values  for  the  ambiguity  vectors  are: 


*1 = 

-.10588E+0£T 

-.19412E+00 

.  .82353 E+00. 

+  ) 

"  .37647E+O0- 

•23529E-01 

-.70588  E+00. 

^2  = 

'  .92496E+0ff 

.28630 E+00 
-.27732  E+00. 

♦  j 

’  .  62643  E+0ff 

.47961 E+00 

-.29364  E+00. 

<7. 16a) 


<7. 16b) 


The  measurements  corresponding  to  actual  parameter  values  are: 


„  [  . 66358 E+Otf 

y  =  .80518E-01  +  j 


M  r-.1A926E+0(n  r-.38510E+0(T 

y  <S2}  =  .38978 E+00  +  j  .22644E+00 


-. 63077 E+Uu | 
.21054E+00j 


(7.17a) 


(7.17b) 


The  Jacobian  for  this  problem  has  dimension  20  x  17  and  therefore  the 
solution  algorithm  requires  the  use  of  the  modified  Newton-Raphson  iteration 
discussed  at  the  beginning  of  this  chapter.  This  is  accomplished  by  the 
IMSLC16]  routine  LLSQF.  The  diagnosis  program  required  five  iterations  and 
1.73  seconds  on  the  VAX  11/780  to  find  r*  to  .1%  accuracy. 

Example  7.18:  Suppose  the  circuit  of  Example  7.9  is  cascaded  to  form 
the  multistage  filter  shown  in  Figure  12.  Diagnosis  of  this  circuit  re¬ 
quired  test  points  between  stages  of  the  amplifiers.  The  reason  is  that  the 
individual  stages  are  isolated  so  that  using  output  measurements  only,  the 
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;■ 


changes  in  a  component  in  one  stage  are  indistinguishable  from  similar 

changes  in  the  corresponding  component  in  another  stage.  The  CCN  equations, 
~R 

1*2^,  and  V  for  the  circuit  of  Figure  12  are  given  below.  They  present  an 
excellent  example  of  the  sparseness  typical  of  the  matrices  of  the  CCM. 
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(7.19) 
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(7.20) 
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u(s) 


u(s) 


MICROCOPY  RESOLUTION  TEST  CHART 
NATIONAL  BUREAU  OF  STANDARoVl963-A 
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V 


1  0  0  0  0  0  0  0 

0  1  0  0  0  0  0  0 

0  -1  000000 
0  0  1  0  0  0  0  0 

00000000 
0  0  0  1  0  0  0  0 

0  0  0  0  1  0  0  0 

0000  -1  000 
0  0  0  0  0  1  0  0 

00000000 
0  0  0  0  0  0  1  0 

0  0  0  0  0  0  0  1 

0000000  -1 
00000000 
00000000 
00000000 
00000000 
00000000 


For  this  example  test  frequencies  are 
solution  parameter  vectors  are: 


0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
1  0 
0  0 
0  1 
0  0 
0  Q 


*  jl  and  s2  *  j2. 


(7.23) 
The  nominal  and 


rQ  *  <2.  1.5  .5  2.5  1.  2.  1.5  .5  2.5  1.  2.  1.5  .5  2.5  1.  .7  1.  2.)* 


(7.24) 

r*  *  (2.1  1.4  .5  2.5  1.  1.4  1.  .5  2.4  1.  1.8  1.6  .45  2.6  1.  .8  1.2  1.8)* 


(7.25) 


The  nominal  values  for  the  ambiguity  vectors  are: 


-■  i  gate 
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'  .49466E+0tf 

'  .978236-01" 

-.103236+00 

-.40269E-01 

-.58126E-01 

-.58Q63E-01 

.123146-01 

.204976-01 

.47561 E-02 

♦  i 

.13267E-01 

-.51815E-03 

-.249836-02 

-.51469E-04 

-.95655E-03 

-.92643E-04 

-.17218E-02 

(7.27b) 


The  diagnosis  program  determined  all  components  of  r*  to  .02%  accuracy. 
This  required  six  iterations  and  31.2  seconds  on  the  VAX  11/780. 

Example  7.18  illustrates  the  fact  that  it  may  be  extremely  costly  in 
terms  of  test  points  (8  required  by  that  example)  to  determine  the  value  of 
every  individual  component.  In  many  instances  individual  components  will  be 
part  of  a  functional  block  which  is  replaceable  as  a  unit.  In  this  situa¬ 
tion  it  is  merely  necessary  to  determine  if  the  functional  block  is  operat¬ 
ing  properly.  The  next  two  examples  illustrate  this  type  of  tradeoff. 

Example  7.28:  Figure  13  depicts  a  twelve  component  system.  Suppose  y^ 
is  a  functional  output  which  makes  it  a  natural  choice  as  a  test  output  as 
well.  From  the  trial  and  error  procedure  described  at  the  beginning  of  the 
section  one  finds  that  the  identification  of  all  the  parameter  values  re¬ 
quires  four  additional  test  points,  y ^  through  y^.  The  CCM  equations, 
and  V  below  reflect  this  choice  of  test  points: 


Figure  13.  System  for  example  7.28 


96 


Let  the  nominal  and  actual  parameter  vectors,  rQ  and  r*  respectively, 
following: 


r0  *  (1  0  1.25  2  -1  .67  .9  .1  2  1.1  -.1  I)1 


r*  *  (.9  .1  1.11  1.8  -.8  .71  1  0  1.33  1.05  -.05  1.05) 


The  nominal  values  for  the  ambiguity  vectors  are: 


'  .165086+0(7 

"  .53947E-0f 

•60717E+00 

.92653E-01 

-.52401 E-02 

-.10637E+00 

•20590E+00 

♦  j 

-.23701 E+00 

-.37403E+00 

-.12865E+00 

. 23949E+00 

-.30979E-01 

-.53498E+00 

-.30480 E+00 

(7.33) 
be  the 

(7.34) 

t 

(7.35) 

(7.36a) 
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‘  . 28803 E+0tf 

‘  .93115E-01* 

.61570E+00 

.23719E+00 

.42695E-01 

-.11083E+00 

•20456E+00 

♦  j 

-.29017E+00 

-.19671 E+00 

-.12473E+00 

•98183E-01 

-.1441 IE-02 

-.30755  E+00. 

-.17381  E+00. 

(7.36b) 


The  measurements  corresponding  to  actual  parameter  values  are: 


H,  * 

y  (s^ 


y"<V 


-.18201 E+0ff 

-.76336E+tKf 

.34693 E+00 

.921 98E-01 

.52903E+00 

♦  ) 

-.39824E+00 

. 78978 E-01 

-.27258E+00 

.  .30958E+00. 

-.47566E+00 

'  .45295E-01" 

-,67300E+0tf 

.56059E+00 

.15584 E+00 

•10124E+01 

♦  j 

-.11559E+00 

.13830 E+00 

-.26937E+00 

.  .32849E+00. 

-.53376E+00. 

(7.37a) 


(7.37b) 


The  diagnosis  program  required  only  three  iterations  and  7.2  seconds  on  the 
VAX  11/780  to  determine  the  elements  of  r*  to  three  significant  digits.  The 
cost  in  terms  of  the  number  of  test  points  is  relatively  high  (five  test 
outputs).  In  the  next  example  the  requirement  for  complete  identification  is 
relaxed  causing  a  dramatic  reduction  in  test  points. 

Example  7.38:  Now  suppose  the  components  of  Example  7.28  permit  func¬ 
tional  grouping  as  indicated  by  the  dashed  lines  in  figure  13.  If  the  re¬ 
quirement  to  isolate  faults  to  a  specific  component  is  relaxed  to  the  iden¬ 
tification  of  a  faulty  group  it  becomes  possible  to  eliminate  test  points  by 
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transforming  each  group  of  three  components  into  individual  components.  For 
example  the  components,  Z^,  z^,  and  Z^,  combine  to  form  the  single  com¬ 
ponent,  Z^  where: 


Zt  (K,  )  = 


"1 

tTs+T 


(7.39) 


and 


K 


1 


r1 

r1+r2 


(7.40a) 


T1 


1 

r1r3+r2r3 


(7.40b) 


If  the  remaining  components  have  similar  transformations  the  twelve  parame¬ 
ter  system  of  figure  13  becomes  the  eight  parameter  (4  component)  system  of 
figure  14.  For  the  system  of  figure  14  the  system  output  y 1  at  four  dis¬ 
tinct  sinusoidal  test  frequencies  is  sufficient  to  identify  the  parameters, 
K.  and  t^,  i=1,..,4  and  hence  isolate  the  fault  to  a  component  group.  The 
CCM  equations  for  the  system  are: 


4 


'  0  0  0-1  " 

V 

T 

0  0  11 

b2 

1 

0  10  0 

b3 

+ 

0 

- 1 

O  O 

O  O 

r-  r* 

r>  t- 

b4 

0 

(7.41) 
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-*  7* 

VI 
♦ 

o 

o 

o 

_ 1 

V 

K2 

V 

b2 

0  — frr  0  0 

TgS+l 

a2 

S 

b2 

K3 

a3 

b4 

0  0  TjS*1  0 

a4 

0  0  0  r^T 

The  parameter  vector  to  be  identified  is: 

The  test  frequencies  are  -  jl,  s^  *  j3,  s^  =  j6  and  =  jlO.  The 
nal  and  actual  parameter  vectors  are: 

rQ  »  (1  2  .9  1.1  .8  1.5  .5  D* 


r*  -  (.9  1.8  1  1.05  .9  1.4  .75  .95)* 


The  nominal  values  for  the  ambiguity  vectors  are: 


’  .  4789805  E+0(f 

’-.6195888E-01" 

“1  “ 

-. 3740277 E+00 
-.5349754  E+00. 

+  j 

-.1286507E+0Q 
- .  3048038E+QG. 

'  .9605259E-0f 

’  .3421843E-0? 

Oj,  * 

-.1289079E+00 

♦  j 

-.8045628E-01 

-.1911152E+00 

-.9222372E-01_ 

(7.42) 


(7.43) 

nomi- 

(7.44) 

(7.45) 

(7.46a) 

(7.46b) 
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'  .2229586E-01" 

’  .7406399E-0? 

—3  a 

-.5140581 E-01 
-.666871 2 E-01. 

♦  j 

-. 20906 79E-Q1 
-.1981 173  E-01. 

‘  . 7343345 E-0? 

'  .6553531 E-Of 

«4  = 

-.21 64005 E-01 
-.2659364  E-01. 

+  j 

-.5590791 E-02 
-.501 0533  E-02. 

(7.46c) 


<7.46d) 


The  measurements  corresponding  to  actual  parameter  values  are: 

yM(s1)  =  -.18201 30E-00  +  j .763361 5e+00  (7.47a) 


y"(s2)  =  .7491249E-01  -  j .5557519E+00  (7.47b) 


yM(s3)  =  .4131773E-01  -  j .3430093 E+00  (7.47c) 

yM(s4)  *  .1806347E-01  -  j .2193868E+00  (7.47d) 

The  program  to  solve  this  example  was  executed  on  the  VAX  11/780  and  re¬ 
quired  9  iterations  to  determine  r*  to  approximately  IX.  The  larger  number 
of  iterations  necessary  for  this  example  versus  the  previous  is  most  likely 
a  result  of  (1)  the  greater  nonlinearity  in  the  tableau  equations  due  to  the 
form  of  the  component  (equation  7.39)  and  (2)  from  the  drastic  reduction  in 
test  points  making  the  extraction  of  parameters  more  difficult. 
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8.  Concluding  Remarks 

This  completes  the  presentation  of  the  preliminary  development  of  the 
tableau  approach  which,  as  pointed  out  earlier,  offers  superior  convergence 
and  conditioning  properties  over  the  equivalent  composite  formulation.  The 
five  examples  demonstrate  the  successful  application  of  the  results  intro-* 
duced  thus  far,  including  the  determination  of  diagnosabi lity  of  a 
circuit/system,  the  calculation  of  the  number  of  required  test  frequencies 
and  the  actual  solution  of  the  fault  diagnosis  equations.  Additional  advan¬ 
tage  for  the  tableau  formulation  remains  to  be  gleaned  from  the  knowledge  of 
the  quadratic  form  of  its  fault  diagnosis  equations.  Details  of  the  exploi¬ 
tation  of  this  information  as  well  as  a  demonstration  of  the  resulting  im¬ 
provement  in  convergence  follow  in  the  next  chapter. 
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CHAPTER  5 

EXPLOITATION  OF  THE  QUADRATIC  FORM 


1 .  Introduction 

This  chapter  introduces  some  improvements  to  the  Newton  algorithm  for 
solving  the  Tableau  Fault  Diagnosis  Equations  (equation  2.2  of  Chapter  A) 
which  apply  when  these  equations  are  quadratic.  Recall  that  the  Tableau 
Fault  Diagnosis  Equations  have  a  quadratic  form  whenever  Z..(s,rO  =  r.Z^(s). 
In  the  case  of  electrical  networks  this  form  of  component  transfer  function 
suffices  for  modelling  resistors,  inductors,  capacitors,  and  dependent 
sources.  This  means  that  the  quadratic  form  of  the  Tableau  Fault  Diagnosis 
Equations  is  valid  for  any  circuit  modelled  by  an  interconnection  of  the 
above  components.  This  is  an  extremely  large  category  of  possible  circuits. 

The  significance  of  this  special  form  is  that  it  represents  specific 
information  about  Tableau  Fault  Diagnosis  Equations  which,  under  appropriate 
circumstances,  is  useful  in  determining  their  solution.  This  chapter 
addresses  two  ways  of  exploiting  this  information.  First  the  quadratic  form 
of  the  Tableau  Fault  Diagnosis  Equations  allows  an  exact  characterization  of 
their  behavior  along  a  search  direction.  In  particular  this 
characterization  provides  a  simple  means  of  precisely  establishing  the  point 
along  the  search  direction  at  which  the  norm  of  the  nonlinear  function, 
F(x),  is  a  minimum.  Use  of  this  "minimum  point"  as  the  next  estimate  of  the 
solution  is  designed  to  guarantee  and  streamline  the  convergence  of  the 
Newton-Raphson  iteration.  Second,  this  special  form  may  permit  the 


determination  of  a  better  search  direction  by  allowing  a  two  term 
approximation  of  the  Hessian  of  the  function,  HFlx)!!^. 

2.  Multidimensional  Quadratic  Functions 

Before  proceeding  with  a  simple  example  of  a  set  of  quadratic 
equations,  it  will  be  useful  to  draw  an  analogy  between  the  problem  of 
solving  F(x)  =  e  and  a  functional  minimization  problem.  Recall  the  square  of 
the  Euclicean  norm,  ||  |||,  is  a  functional  which  has  a  unique  global 
minimum  at  6  so  the  problem  of  finding  x*  such  that  F(x*)  *  0  is  equivalent 
to  finding  x*  such  that  C133 

1 1  F<x*)  1 1|<^|  |  FCx)|  ||  for  all  xeRn  (2.1) 

For  circumstances  which  cause  the  Newton-Raphson  algorithm  to  diverge,  this 
new  perspective  suggests  a  modification  which  results  in  a  nondive’ging 
algorithm  C133.  Namely,  choose  x  such  that 

||F(xk+1)|||  »  ||F(xk+xdk)|||  <  ||F(xk)|||  (2.2) 


|r 

where  d  satisfies: 

Jf(xk)dk  =  -F(xk)  (2.3) 

The  vector  d  is  called  the  "search  direction".  This  modification  is  based 
on  the  "Global  Convergence  Theorem"  which  appears  in  Chapter  2.  The 
following  example  will  illustrate  how  knowledge  of  the  quadratic  form 
assists  in  the  selection  of  x. 


Consider  the  following  problem  in  two  dimensions 
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F(x,y> 


(2.4) 


where  J^jeR*.  The  nonlinearity  of  the  function,  F,  is  the  second  order 
term,  xy.  Such  second  order  terms  likewise  characterize  the  nonlinearity  of 
the  Tableau  Fault  Diagnosis  Equations. 

Now  solve  F(x,y)  =  e  via  the  Newton-Raphson  iteration  method.  The 
Jacobian  of  F(x,y)  is 


Assume  an  initial  guess: 

r.i 

l/'J 


solution. 


i/j 

,  is  the  solut 


v«->  -  [?  ::g 


(2.5) 

The  next  estimate  of  the  desired 


on  of: 


Je<*6#y°> 


(2.6> 


which  is  equivalent  to  solving  the  pair  of  linear  simultaneous  equations: 

x  +  .1y  *  3.3  (2.7a) 


3x  +  ,1y  =  5.3  (2.7b) 

To  develop  an  understanding  of  the  nature  of  the  solution  process 
consider  the  following  physical  interpretation:  Let  the  first  component  of  F 
define  a  surface  as 


z  *  F1  (x,y> 


(2.8) 
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Figure  IS  is  a  graph  of  this  surface 


with 


the 


point 


z°  *  F^(x°,y°)  =  (2.1/3)  annotated.  Consider  next  the  plane  which  is 

tangent  to  the  surface  at  the  point  x  I,  depicted  in  figure  16.  Equation 

2.7a  is  the  equation  of  the  line  in  the  xy  plane  which  is  the  intersection 

of  the  tangent  plane  with  the  plane  z*0.  In  other  words,  equation  2.7a  is 

the  linear  approximation  of  the  zero  set  of  the  nonlinear  function,  F^(x,y) 

determined  as  the  zero  set  of  the  tangent  plane  of  the  surface  z  =  F1(x,y) 

f>l 

at  the  point  .  Of  course  a  similar  interpretation  holds  for  the  second 
o 

Ly  J 

equation.  Thus  the  intersection  of  the  two  tangent  planes  is  the  approximate 
intersection  of  the  two  nonlinear  equations. 


The  solution  to  equation  2.7  is  the  point 


C3-W- 


To  determine  if 


the  algorithm  may  be  diverging,  test  the  condition  in  2.2,  i.e.  determine 
if  1 1 F Cl ,23) | <  1 1 F(2.1 ,3) 1 1 Evaluating  1 1 F 1 1 ^  at  these  points  yields: 


F <2.1 ,3) 1 1 j  =  2*5  1 1 F (1,23) I Ig  *  968 


(2.9) 


This  indicates  that  the  Newton-Raphson  iteration  is  diverging.  For  this 
simple  example  it  is  very  possible  that  subsequent  iterations  may  recover 
from  this  divergent  step  but  in  general,  adjustment  of  the  algorithm  to 
insure  that  the  condition  in  2.2  holds  will  be  necessary. 

This  adjustment  consists  of  finding  a  value  for  X  such  that 

||F(xk+Xdk)| \\  <  ||F(xk)| |*  (2.10) 

For  this  example  the  search  direction,  d°  =  fc]  -  *  [Ifl  .  The 

desired  value  for  the  next  estimate  of  the  solution  is  now 


F*wW. 
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rfl 

’  0 

d° 

X 

2 

X 

X  1 

X 

l/’J 

A 

▼  A 

[4 

(2.11) 


Notice  that  when  x  =  0  the  next  estimate  equals  the  previous  and  when  x  *  1 


the  next  estimate  Is  the  Newton-Raphson  estimate.  The  set  of  possible 

:j\ 

1  .  corresponds  to  a  line  parallel  to  the  search  direction  and 
1  -ol 

.  Figure  17  Illustrates  this  set  of  points  for  the 


values  for 

l/'J 

passing  through 


current  example.  In  the  Navld  and  Willson  application  of  the  "Global 
Convergence  Theorem"  x  has  an  Initial  value  of  1  (the  Newton-Raphson 
Iterate).  If  the  condition  In  2.2  does  not  hold  then  x  decreases  by  half 
until  It  does  C13D.  Points  corresponding  to  x  *  1,.5  and  .25  appear  In 
figure  17.  Table  2  displays  the  results  of  the  computation  of  ||F|||  for 

tk.  ianro  ?  ) 


for  j=0,1,2 

and  3. 

X 

imi| 

1 

968 

.5 

58.93 

.25 

4.363 

.125 

1.910 

Table  2 

i2 


Evaluation  of  | | F | | ^  at  selected  points. 


The  data  In  table  2  Indicates  that  x  =  .125  will  satisfy  the  condition  In 
2.2.  Although  this  value  for  x  satisfies  the  condition.  It  Is  not  optimal  In 
the  sense  that  IlFll^  Is  not  minimized  along  the  search  direction. 
Knowledge  of  the  quadratic  nature  of  the  Tableau  Fault  Diagnosis  Equations 
and  hence  the  fourth  order  behavior  of  ||F|||  offers  an  easily  Implemented 
scheme  for  choosing  this  optimal  value  as  seen  In  the  next  section. 


Figure  17.  Set  of  Points  Along  the  Search  Direction. 


Ill 


3.  An  Algorithm  for  Setting  \ 

To  illustrate  an  optimal  scheme  for  choosing  a,  return  to  the  example 
in  the  previous  section.  Evaluate  F(x^,y^)  as  a  function  of  A  by 
substituting  equation  2.11  into  the  expression  for  F(x,y)  given  in  equation 
2.4.  This  produces: 


F (2.1~1 .1 A,3+2QA)  *  |”22x  +*9X"'9  1  (3.D 

L-22X  -1.3X+1.3J 

Note  that  F  is  a  vector  whose  entries  are  quadratic  functions  of  X.  Now 
find  A  such  that  ||F(2.1-1.1A,3.+20.X)||2  is  a  minimum.  Computing  the 
square  of  the  norm  of  F  yields  the  following  fourth  degree  polynomial  in  A: 

||F<x1,y1)|||  *  <-22A2+.9X-.9>2  ♦  C-22a2-1 .3A+1 .3>2 

*  968A4+17.6A3-15.1A2-5X*2.5  (3.2) 


Since  equation  3.2  is  a  fourth  degree  polynomial  the  minimum  point  is  the 
result  of  solving: 

d||F(x1,y1>||?  3  2 

- - -  =  3872A  t52.8a  -30. 2A-5  *  0  (3.3) 

Equation  3.3  has  a  single  real  root  at  A  *  .12731.  Thus  the  point 

is  the  point  along  the  search  direction  for  which  IlFll^  is  a  minimum. 
Evaluating  F  at  this  point  gives 
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1 1 F<1 .96,5.5462) |||  *  1.909 


(3.5) 


This  process  has  produced  a  new  estimate  for  the  solution  of  F(x,y)  *  0  void 
of  the  divergence  problem. 

The  procedure  for  finding  the  optimum  value  for  x  worked  very  well  for 
the  simple  example.  The  next  step  in  this  development  is  to  show  explicitly 
that  it  applies  to  the  Tableau  Fault  Diagnosis  Equations.  Recall  the  fault 
diagnosis  equations  presented  in  equation  3.1  of  chapter  3,  namely: 


F  (x)  = 


‘f1<r)g1(a1)-0i* 


,f  <r)g  (a  )-0  , 

L  q  n  q| 


*  e 


(3.6) 


where 


<i) 

(ii) 


(iii) 


x  s  co l  (j^| /<i./ . .  > , r )  $ 
f • (r)  =  CZ(s.,r) |-V3  ; 

P-iVW’i'l 


9,(5,) 


a. 

—1 


;  and 


(iv) 


®i  “  bo<si)* 


Furthermore  recall  the  assumption  that  Z^(s,r^)  =  r^Z^(s). 


To  explicitly  see  that  the  procedure  used  in  the  simple  example  applies 
to  these  equations  as  well,  take  a  closer  look  at  the  individual  entries  of 
the  nonlinear  vector  function,  F.  Let  $„(x)  denote  the  CN(i-1)+j3th 
component  of  F(x),  (i.e.  FN(i-1)+j*x^  where  is1,2,...,q  and  js1,2,...,N. 


Then 
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♦  z 


i'W.i'V-Eyirti'v 


(3.7) 


where 

(i)  (L^VK^  denotes  the  (j,k)  entry  of  the  matrix  product  L^V; 

<ii)  o^k  is  the  k-th  entry  of 
(iii)  denotes  the  <j,k)  entry  of  the  matrix  V;  and 

(iv)  a  . (s.)  and  b  ,(s4)  denote  the  j-th  entry  of  a  (s.)  and  b  (s.) 

respectively. 

The  nonlinear  terms  are  those  with  the  products:  r.o..  where  i=1,2,...,q, 

J  I K 

j*1,2,...,n,  and  k=1,2,...,p.  Clearly  then  the  nonlinearity  in  the  Tableau 
Fault  Diagnosis  Equations  consists  of  sums  of  pairwise  products  of  r.  with 
a.k,  strictly  quadratic. 

Now  apply  the  modified  Newton-Raphson  iteration  to  the  Tableau  Fault 
Diagnosis  Equations  in  the  same  manner  as  in  the  earlier  simple  illustrative 

example.  As  before  the  objective  is  to  select  the  value  of  x  such  that 

k  k  2  k+1  k  k 

| | Fix  +xd  > | | 2  is  a  minimum  and  then  compute  x  ■  x  +Xd  .  From  equation 

k  k 

3.7  it  is  apparent  that  with  x  and  d  fixed  (known) 
k  k 

(i)  ♦^(x  +xd  )  is  a  quadratic  function  of  X,  and 

lii)  ||F(xk+Xdlc)|  | g  is  a  fourth  degree  polynomial  in  X. 

Furthermore  the  above  is  true  regardless  of  the  size  of  the  Tableau  Fault 
Diagnosis  Equations  as  long  as  Z^ls,^)  »  rjZjls). 

This  means  that  the  following  algorithm  will  assign  the  optimum  value 

to  x: 


| 

I 
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Step  1-  Evaluate  1 1  F(xk*xdS  1 at  five  values  of  X. 

Remark:  Because  of  the  form  of  the  Tableau  Fault  Diagnosis  Equations,  this 
requires  matrix  addition  and  multiplication  only.  No  matrix  inverses  are 
necessary. 

Step  2-  Generate  a  fourth  order  interpolant,  h(x),  for  the  data  from 
Step  1. 

Remark:  Although  it  is  possible  to  generate  the  coefficients  of 

k  k  ?  k  k 

| I F(x  +xd  )| |-  by  substituting  x  +Xd  into  the  expression  for  F  as  in  the 
example  above  it  is  far  simpler  in  a  computer  implementation  to  use  steps  1 
and  2  for  the  following  reasons.  First,  any  computer  implementation  of  the 
Newton-Raphson  algorithm  (e.g.  the  program  used  to  solve  the  examples  in 
chapter  5)  must  already  have  a  routine  for  evaluating  the  value  of  F. 
Second,  since  the  | | F | | 2  is  a  fourth  degree  polynomial  in  a  then 

h(A)  =  | | F<xk+Adk> | | |  (3.9) 


for  all  x  (i.e.  the  interpolant  is  exact). 
Step  3-  Find  the  roots  of  the  equation: 


dh(A)  _ 
dx  " 


(3.10) 


Remark:  The  solution  to  equation  3.10  is  easily  computed  since  is  a 

third  degree  polynomial. 

Step  4-  Set  x^  equal  that  root  of  equation  3.10  for  which  3.9  is 
minimum. 

The  result  of  this  procedure  is  a  value  for  x=xk  which  minimizes 
||F(xk-xdk)| ||  along  the  search  direction,  dk.  This  is  the  optimum  choice 
for  x  made  possible  by  the  quadratic  nature  of  the  Tableau  Fault  Diagnosis 
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Equations. 

Before  proceeding  to  an  example  some  details  of  the  implementation  of 
steps  1  and  2  deserve  mention.  If  infinite  precision  were  available  on  the 
computer  any  five  distinct  points  along  the  search  direction  would  suffice 
to  uniquely  determine  the  desired  interpolant,  h(x>.  With  limited  precision 
some  choices  will  produce  more  accurate  results  than  others.  Since  the 
solution  to  equation  3.10  will  normally  lie  in  the  interval  CO, ID  a 
reasonable  choice  of  five  interpolation  points  is  x  =  -.5,0. ,.5,1 .,  and  1.5. 

The  polynomial  h(X)  has  the  form: 

h(x)  =  y4X4  +  YjX3  +  y2x2  +  ♦  Y0  (3.11) 

Evaluation  of  h(X)  =  | | F(xk-Xdlc) | |^  at  the  five  points  mentioned  above 
results  in  the  following  matrix  equation  whose  solution  is  the  coefficients 
Y.j,  i=0,1,..,4,  which  uniquely  specify  h(X). 


[5.0625  3.375  2.25  1.5  1.1 

[h<1 .5)"] 

1.  1.  1.  1.  1. 

Y3 

h(1.0) 

0.0625  0.125  0.25  0.5  1. 

y7 

= 

h(0.5> 

0.  0.  0.  0.  1. 

c 

h(0.0) 

0.0625  -0.125  0.25  -0.5  1 .. 

Ti 

h(-.5) 

;o_ 

The  matrix  in  equation  3.12  is  a  5  x  5  Vandermonde  matrixCKD  which  is 
nonsingular.  Since  this  matrix  is  always  the  same  as  long  as  the 
interpolation  points  are  fixed  the  solution  algorithm  computes  its  inverse 
once  and  then  uses  the  inverse  at  each  iteration  to  determine  the 
coefficients  of  h(X). 
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4.  Quadratic  Fora  Examples 

This  section  will  illustrate  the  modified  Newton-Raphson  algorithm 
(selection  of  x)  detailed  in  the  previous  section  via  an  example  circuit 
(Figure  9)  presented  in  chapter  4. 

Example  4.1:  The  first  example  is  example  7.1  from  chapter  4.  All  the 
data  used  in  this  example  is  exactly  the  same  as  in  example  7.1  except  that 
this  time  the  solution  to  the  Tableau  Fault  Diagnosis  Equations  proceeds 
according  to  the  modified  Newton-Raphson  Algorithm.  Recall  that  for  this 
example  M=N=12,  p=9  and  q=2.  The  test  frequencies  are  s^=j1.4  and  s^=j2.0. 
Determination  of  the  solution  point  with  an  accuracy  of  .1 1 X  required  8 
iterations  and  23.3  seconds  on  the  VAX-1 1/780.  The  solution  routine  is  a 
FORTRAN  language  program  using  single  precision  arithmetic  (approximately 
six  decimal  places  on  the  VAX).  This  is  a  substantial  improvement  over  the 
Newton-Raphson  iteration  (xs1)  described  in  chapter  4,  example  7.1,  which 
required  12  iterations  and  31.5  seconds  on  the  VAX-11/780  to  produce  the 
same  results. 

Example  4.2;  This  example  also  uses  the  same  circuit  (Figure  9)  as  the 
previous  example  with  the  exception  of  the  actual  parameter  vector  and 
corresponding  measurement  outputs.  For  this  example  the  actual  parameter 
vector  is: 

r  *  (2  1.4  .9  .7  1.5  .8  1.5  1.4  2  .9  1.2  1.6)*  (4.23) 


and  the  test  outputs  corresponding  to  this  parameter  vector  are 


-.14522E-0f 

\29083E-0f 

.25990E+00 

♦  j 

. 54699E+00 

_  .65868E-01. 

.13866E+0Q 

(4.4a) 


for  input  u(J1.4)  *  1,  and 


r.  1 2406E-0ifl  r.471 75E-01" 

yM(j2.0)  =  .43757E+00  ♦  j  .66102E+00  <4.4b) 

94981  E-OlJ  1 8066E+0Q 

for  Input  u(j2.0)  =  1. 

The  variations  from  unity  for  each  of  the  parameter  values  are  sufficiently 
large  that  the  usual  Newton-Raphson  iteration  diverges.  Consequently  this 
problem  solution  required  the  modified  algorithm  discussed  in  section  3. 
Determination  of  the  solution  point  with  an  accuracy  of  ZX  required  6 
iterations  and  17.0  seconds  on  the  VAX-11/780. 

The  examples  illustrate  that  the  modified  algorithm  will  converge  under 
circumstances  in  which  the  Newton-Raphson  algorithm  will  not.  Furthermore  it 
is  possible  for  the  modified  algorithm  to  result  in  convergence  in  fewer 
iterations  due  to  the  optimum  choice  of  X.  Of  course  these  conclusions  are 
consistent  with  the  "Global  Convergence  Theorem"  of  Navid  and  Willson  C133. 

5.  Modification  of  the  Search  Direction 

The  algorithm  for  solving  the  quadratic  form  of  the  Tableau  Fault 
Diagnosis  Equations  uses  this  quadratic  information  at  each  iteration  to 
determine  the  optimum  point  along  the  search  direction  for  the  next  estimate 
of  the  solution.  The  search  direction  used  however  is  exactly  the  same  as 
the  Newton-Raphson  algorithm.  In  both  cases  the  search  direction  for  the 

|r 

k-th  Iteration,  d  ,  satisfies: 

JF<xk)dk  =  -F(xk>  (5.1) 

The  purpose  of  this  section  is  to  Investigate  the  use  of  quadratic 
information  to  determine  of  a  new  possibly  better  search  direction. 
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In  section  2  the  problem  of  solving  the  nonlinear  vector  problem, 
FCx)  =  e,  appears  recast  as  a  the  minimization  of  a  scalar  function 
1 | F (x) | | ^  where  xeRn.  Let  the  function,  h(x),  be  a  scalar  multiple  of 
| | F(x)| ||,  expressed  in  terms  of  vector  multiplication  as: 

h(x)  =  .5  F(x)tF(x)  (5.2) 


where  F:  DC  Rn->R*  with  m>n  and  t  denotes  the  matrix  transpose.  The  scale 
factor  .5  has  no  affect  on  the  minimization  problem  but  its  inclusion 
eliminates  the  factor  two  from  the  derivatives  of  h(x).  Recall  that  finding 
the  solution  to  F(x)  *  0  is  equivalent  to  finding  the  value  of  x  which 
minimizes  h(x).  Furthermore,  finding  the  minimum  of  h(x)  is  equivalent  to 
solving  the  nonlinear  vector  problem: 

VhCx)  =  e  (5.5) 


where  v  denotes  the  gradient  operator,  i.e. 


v  s 


3x„ 


»x. 


8xn 
L  J 


(5.4) 


Define  the  function  G(x)  as: 


G(x)  »  vh(x) 


(5.5) 
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The  Nevton-Raphson  and  modified  Newt on-Raph son  algorithms  are  valid 
methods  for  finding  the  solution  to  6(x)  =  e.  In  both  cases  the  search 
direction  at  the  k-th  iteration  is  the  solution  to: 

J6<xk)dk  =  -G(xk>  (5.6) 

Expressing  equation  5.6  in  terms  of  the  original  nonlinear  function,  F(x), 
produces  an  interesting  variation  on  the  Newton-Raphson  search  direction.  In 
order  to  see  this  take  a  closer  look  at  the  Jacobian  of  G.  If  the  vector 
function 

G(x)  *  colCG. (x),G0(x),.  .  .,G  (x)3  (5.7) 

■  c  n 

_ ,  fah(x)  ah(x)  8h(x)l 

-  coiL— J 

then  J^(x)  is 
6 
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d2h(x) 

d2h(x) 

32h(x) 

3X13X1 

3x23x1  * 

*  3xn3X, 

32h(x) 

32h(x) 

a2h(x) 

3x.j3X2 

0 

3X23X2  * 

•  • 

•  9 

9 

32h(x) 

•  • 

32h(x) 

0  0 

a2h(x) 

3x.  3x 

1  n 

3x,3x 

2  n 

•  ix  lx 
n  n 

The  v2  operator  denotes  the  Hessian  matrix  (the  matrix  of  second  partial 
derivatives. 

Both  the  gradient,  G(x)  *  vh(x),  and  the  Hessian  matrix,  ?2h(x),  have 
expressions  in  terms  of  F(x)  based  on  the  definition  of  h(x)  given  in 
equation  5.2.  These  are  C223: 


G(x)  »  Vh(x)  * 


(5.9) 


v2h(x)  = 


t  * 

JP(x)  ♦  £  F.(x)H.(x) 
f  k»i  K  * 


(5.10) 


where  Fk<x>  is  the  k-th  entry  of  F(x)  and  Hk(x>  is  the  Hessian  matrix  of 
Fk(x)  i.e. 

Hk(x)  *  vZFk<x)  «.11) 

The  substitution  of  the  identities  of  equations  5.9  and  5.10  into 
equation  5.6  produces  the  following  equation  for  the  search  direction: 
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(5.12) 


Omitting  the  term 


m 


2  F.(x)H.(x>  from  equation  5.12  yields 
k=1  K  * 


(5.13) 


Solving  equation  5.13  is  mathematically  equivalent  to  finding  the  least 
squares  solution  to  equation  5.1  which  produces  the  Newton-Raphson  search 
direction.  Equation  5.12  on  the  other  hand  produces  a  search  direction  based 
on  a  complete  expression  of  the  Hessian  matrix  of  h(x).  An  algorithm  based 
on  the  use  of  equation  5.12  might  therefore  produce  an  improved  solution 
scheme  for  solving  F(x)  *  e.  As  will  be  seen  from  the  examples  of  the  next 
section,  the  procedure  does  not  in  fact  produce  an  improved  search 
direction. 

m 

In  general  information  to  permit  the  computation  of  £  F. (x)H  (x)  (i.e. 

k=1  *  k 

the  Hfc)  is  not  available  and  a  great  variety  of  schemes  exist  which  attempt 
to  approximate  it  in  some  manner  C22D.  In  the  case  of  the  quadratic  form  of 
the  Tableau  Fault  Diagnosis  Equations  the  matrices  H^lx)  =  v^f^(x), 
k=1,2,...,m  are  constant.  Furthermore  they  are  sparse  which  permit 
efficient  storage  for  use  In  computing 

^  t  m 

(x)  J.(x)  ♦  F.(x)H 

J  p  k*1  * 

at  each  iteration.  The  sparseness  is  due  to  the  fact  that  the  only  nonlinear 
terms  in  F^  are  the  products  of  a  single  component  parameter,  and  the 
components  of  a  single  ambiguity  vector,  a.,  i»1,2,...,q  and  j«1,2,...M. 
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Specifically,  each  is  a  <pq+N) (mu(pq+N)  symmetric  Matrix  which  has 

all  zero  entries  except  for  the  (pq+i)-th  row  and  column.  If  r  _  . 

Pd+)*i 

denotes  the  non-zero  column  of  HucMHj  then  the  non"*ero  row  i*  Fpq+j  •{* 
r _ . ■  .  has  the  form: 


Pd+j/i 


(5.14) 


where 


(i) 

equals  the  p  dimensional  zero 

vector 

if  k#i  and 

Ak  *  Zj (Sj)  (L^V)*,  if  K  *  i; 

(ii) 

(L^V)^  is  the  transpose  of  the  j-th  row 

•  4 
> 

-J 

o 

and 

(iii) 

eN  is  the  N  dimensional  zero  vector. 

The  results  of  the  implementation  of  this  modified  search  direction 
method  are  discussed  in  the  next  section. 


6.  Examples  Using  Modified  Search  Direction 

As  with  the  other  solution  algorithms  discussed,  the  implementation  of 
the  solution  algorithm  utilizing  the  modification  of  the  search  direction  is 
a  FORTRAN  program  compiled  and  executed  on  the  VAX  11/780.  The  program  also 
employs  the  interpolation  procedure  to  determine  the  optimum  point  along  the 
search  direction. 

Example  6.1_:  Since  example  4.1  has  already  served  as  a  basis  for 
comparison  between  the  Newton-Raphson  algorithm  and  the  modified  Newton- 
Raphson  algorithm  discussed  in  section  3  it  is  appropriate  to  attempt  to 
solve  this  example  utilizing  the  modified  search  direction.  Unfortunately 


5 

n 


A 
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the  results  were  discouraging.  After  250  iterations/  the  program  utilizing 
the  modified  search  direction  failed  to  converge  to  a  solution  point.  The 
following  conjecture  may  serve  as  an  explanation  of  this  unsatisfactory 
performance.  Recall  that  the  modified  search  direction  algorithm  is 
actually  a  procedure  to  solve  vh(x)  *  e.  This  condition  is  satisfied  not 
only  at  the  desired  solution  point  but  at  all  local  maxima  and  minima  and  at 
saddle  points.  In  fact  it  is  possible  that  some  of  the  points  satisfying 
Vh(x)  =  e  are  not  even  isolated.  By  construction  the  modified  search 
direction  is  influenced  by  all  such  points  on  the  surface  described  by 


The  existence  of  non* Isolated  solutions  or  even  clustered  isolated 


solutions  to  vh(x)  =  e  would  explain  the  program's  inability  to  converge  to 
a  solution  point,  even  an  incorrect  one. 

There  is  perhaps  a  very  limited  set  of  situations  in  which  use  of  the 
modified  search  direction  is  appropriate  as  the  next  example  illustrates. 

Example  6.2:  Like  examples  4.1,  4.2,  and  6.1,  this  example  employees 
the  circuit  of  example  7.1  of  chapter  4.  All  of  the  necessary  data  appears 
there  except  for  the  following  changes: 

(1)  The  number  of  test  f requencies,q,  is  three,  and  these  are 
s1  *  j .4,  s1  =  jl .2,  and  s^  *  j2.; 

Remark:  This  example  uses  one  more  test  frequency  than  actually  required. 
The  additional  measurement  data  serves  to  improve  the  condition  of  the 
problem. 

(1i)  The  actual  parameter  vector  Is: 

r*  *  cot (2.1  1.5  .4  1.2  2  .7  .8  1.1  .95  1.2  .95  1.3)  (6.3) 


(ill)  The  nominal  values  for  the  ambiguity  vectors  are: 


If 
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‘  .26648E+O0 

-.12054E+0CT 

.37331 E+00 

-.11670E+00 

.35839E-01 

.43052E-01 

. 28384 E +00 

-.43458E-01 

-.17221 E-01 

♦  ) 

.14336E-01 

.54351 E-01 

. 1 7588E-01 

.35839E-01 

.43052E-01 

.31641 E-01 

.97458E-01 

.92483E-01 

.14854E+00 

atSj) 


«(Sj) 


r  .12731  E+Otfl 


-.38703E+0dl 


.25786E+00 

.14630E-02 

.86319E-01 

.34182E-01 

.18968E+00 

-.49547E-01 

-.41018E-01 

+  j 

.10358E+00 

.72697E-01 

-.80235E-01 

.86319E-01 

.34182E-01 

.15326E+00 

.11745E+00 

_  .23754E+0Q 

_  . 20604 E+00_ 

'-.86316E-0f 

-.48784E+0ff 

. 29425 E+00 

•80346E-01 

•10092E+00 

.10641 E-01 

.17635E+00 

-.11167E-01 

-.21282E-01 

+  ) 

.20184E+00 

.35081 E-01 

-.13916E+00 

.10Q92E+00 

.10641E-01 

.19156E+00 

.97295E-01 

.31972E+00 

.21512E+00 

(6.4b) 


(6.4c) 


and 

(iv)  The  outputs  corresponding  to  the  actual  parameter  vector  are: 


„  [-.26178E-0ifl  {-.31569E-07 

y  (s,)  =  .31077E-01  ♦  j  .18086E+00  (6.5a) 

1  L  .2111 8E— OlJ  L  .57588E-01 
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-  f-.20198E-0f|  f  .17205E-01" 

y"<s.)  *  .20012E+00  +  j  .45228E+00  (6.5b) 

*  .64193E-01  L  -12685E+00 


„  P.78967E-021  T  .48328E-0f| 

y"<sT>  »  .41679E+00  +  j  .59643E+00  (6.5c) 

3  .11324E+00J  L  .17160E+00J 

The  solution  of  the  resulting  fault  diagnosis  equations  was  coaputed 
twice.  The  first  solution  program  used  the  modified  Newton-Raphson  algorithm 
developed  in  section  3.  This  program  required  11  iterations  and  66.13 
seconds  on  the  VAX  11/780  to  compute  the  solution  to  .07X  accuracy.  The 

second  solution  program  was  identical  to  the  first  except  that  the  first 

iteration  used  the  modified  search  direction  discussed  in  section  5.  This 
solution,  computed  to  .07*  accuracy,  required  8  iterations  and  58.07  seconds 
on  the  VAX  11/780.  In  this  case  a  modest  improvement  resulted  from  the  use 
of  the  modification  to  the  Newton-Raphson  search  direction. .  Experience 
indicates  that  such  improvement  does  not  occur  in  general. 

7.  Summary 

Clearly  a  notable  distinction  of  the  use  of  the  Tableau  Fault  Diagnosis 
Equations  is  the  ability  to  exploit  their  special  structure  in  the  solution 
process.  The  most  significant  observation  about  this  structure  is  that  the 
equations  are  quadratic  when  Z^(s,r^)  *  r^Z.(s).  This  means  that  it  is 
possible  to  completely  characterize  the  behavior  of  ||F|||  in  a  search 
direction  with  a  simple  Interpolation  scheme.  The  interpolation  scheme 
combined  with  the  Newton-Raphson  search  direction  produces  a  solution 
algorithm  which  converges  under  conditions  in  which  the  Newton-Raphson 
algorithm  does  not  or  is  faster  under  conditions  when  it  does. 


Although  the  attempt  to  develop  an  Improved  search  direction  based  on 
the  use  of  a  coaplete  expression  for  the  Hessian  of  1 1 F 1 1 2  was  unsuccessful 
it  does  provide  additional  infornatlon  into  the  nature  of  the  Tableau  Fault 
Diagnosis  Equations.  In  general  the  preserve  of  extrena  and  other  points 
for  which  vh(x>  *  e  precludes  the  use  of  this  aethod.  In  fact  this  should 
serve  as  a  caution  to  the  application  of  this  Method  to  the  solution  of  any 
set  of  similtaneous  nonlinear  equations. 
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